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impossible^investigations  of  the  relationships  of  the  dynamic 
three-dimensional  structural  geometry  to  the  physical  and  biochemical 
functions  of  both  moving  and  stationary  organ  systems. r 
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PROTECTION  OF  THE  CARDIOPULMONARY  SYSTEMS  AGAINST  THE 
INJURIOUS  EFFECTS  OF  ACCELERATION 


Evaluation  of  effects  of  time  history  of  gravitational- 
inertial  force  environment  on  regional  dynamic  distribution 
of  pulmonary  ventilation,  strains,  edema  and  blood  volume. 

Introduction 


It  is  of  fundamental  importance  in  studies  of  cardiovascular  and 
pulmonary  physiology  to  measure  and  analyze  the  dynamic  changes  in 
regional  geometry  and  perfusion  in  these  organ  systems.  However, 
the  function  and  performance  of  the  heart,  lungs,  and  circulation 
are  inseparably  linked  to  their  motion  in  three  dimensions.  Further- 
more, since  the  function  and  interaction  of  the  heart,  lungs  and 
circulation  are  readily  affected  by  physical  interventions  into 
their  working  environment,  and  since  these  interventions  may  result 
in  functional  compensations  which  obscure  the  normal  functional 
status  or  pathological  condition  of  the  organ  or  organ  system 
being  studied,  it  is  desirable  to  devise  measurement  techniques 
whose  nature  involves  the  minimum  possible  degree  of  physiological 
disturbance. 

The  transmission  of  x-rays  through  the  chest  produces  images 
(shadowgraphs)  with  high  temporal  and  spatial  resolution  without 
subjective  sensations,  and  is  non-invasive  in  that  the  x-rays  do 
not  directly  affect  the  function  of  the  cardiovascular  and  pulmonary 
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systems  at  the  dose  levels  required  to  produce  useful  images. 

Since  x-rays  readily  penetrate  and  are  absorbed  to  some  degree 
by  all  bodily  tissues,  they  can  be  used  to  produce  radiographic 
images  of  every  body  organ  ranging  in  density  from  bone  to  lung. 

Television  techniques  provide  the  capability  of  recording  x-ray 
images  at  high  temporal  resolution  in  electronic  form,  which 
greatly  facilitates  their  transmission,  display  and  more 
importantly,  for  quantitative  dynamic  studies,  their  very  rapid 
and  direct  conversion  into  digital  form.  High-speed,  high- 
fidelity  digital  conversion  of  the  images  allows  the  tremendous 
capability  of  digital  computers  for  handling  the  high  data  rates, 
data  volumes  and  numbers  of  calculations  required  for  dynamic 
analyses  of  image  densities  and  dimensions  to  be  brought  to  bear 
on  important  biological  problems. 

Thus,  the  basic  determinants  required  for  accurate  estimation  of 
the  functional  status  of  the  intact  heart,  lungs,  and  circulation, 
namely,  their  true  dynamic  changes  in  regional  shape  and  dimensions, 
the  temporal  and  regional  distribution  of  blood  flow  to,  within, 
and  from  these  organs,  and  their  internal  and  transmural  pressures, 
can  be  obtained  by  application  of  dynamic  x-ray  video  imaging 
techniques  coupled  with  synchronous  recordings  of  multiple 
associated  physiologic  variables. 

It  is  this  information  that  is  a basic  requirement  for  understanding 
the  factors  which  determine  man's  physiologic  response  and  tolerance 
to  short  and  long  term  changes  in  the  force  environment  such  as 
encountered  in  conventional  aircraft  and  in  manned  space  flight. 
Recognizing  this,  a major  effort  of  our  laboratory  over  the  last 
fifteen  years  has  been  directed  toward  the  development,  imple- 
mentation and  application  of  dynamic  quantitative  imaging  techniques 
to  study  the  heart,  lungs,  and  circulation.  In  particular,  two 
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unique  and  complimentary  computer-based  roentgen-video  techniques 
have  been  developed  (Figure  1)  and  provide  the  heretofore  lacking 
technol<^ic  capabilities  to  study  the  spatial  distribution  of 
lung  parenchymal  stresses  and  strains,  regional  lung  volumes, 
regional  myocardial  perfusion  and  the  true  dynamic  changes  in 
shape  and  dimensions  of  the  intact  heart,  lungs  and  thorax  under 
normal  physiologic  conditions,  namely,  an  intact  thorax  and  normal 
intrathoracic  pressures. 

SPATIAL  AND  TEMPORAL  RELATIONSHIPS  OF 
DYNAMIC  REGIONAL  LUNG  MECHANICS 


Figure  1 The  parenchymal  marker  technique  measures  the 

dynamic  spatial  distribution  of  lung  parenchymal 
strains  and  regional  lung  volumes  while  the  three- 
dimensional  reconstruction  technique  measures  the 
dynamic  changes  in  regional  shape  and  dimensions 
of  the  heart,  lungs,  chest  wall  and  diaphragm. 
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This  final  progress  report  for  Contract  No.  F-44620-71-C-0069, 
covering  the  period  July  1,  1971  to  June  30,  1976,  consists  of 
three  parts  and  an  appendix.  Parts  I and  II  contain  a compre- 
hensive review  and  current  status  of  the  pulmonary  parenchymal 
marker  (Part  I)  and  three-dimensional  reconstruction  (Part  II) 
techniques  and  the  results  derived  therefrom.  Details  of  other 
studies  pursued  during  this  period  have  been  documented  in  published 
papers  and/or  previous  annual  progress  reports.  Part  III  consists 
of  a list  of  publications  from  our  laboratory  which  were  based  on 
work  supported  by  funds  from  this  contract.  The  Appendix  contains 
a manuscript  by  Ralph  Sturm  and  co-workers  in  our  laboratory  which 
describes  in  detail  a new  state-of-the-art  video  imaging  system 
(i.e.,  SSDSR)  which  was  designed  by  them  and  was  used  for  the 
three-dimensional  reconstruction  studies  described  in  Part  II. 

PART  I : Spatial  (i.e.,  Three-Dimensional)  Distribution  of  Regional 

Lung  Strains  and  Volumes  in  the  Intact  Canine  Lung  From 
Measurements  of  the  Dynamic  Displacements  of  Roentgen- 
opaque Parenchymal  Markers  in  Anesthetized  and  Awake  Dogs 

Under  dynamic  conditions,  the  mechanical  properties  of  the  lungs 
and  thorax,  including  the  diaphragm,  play  an  important  role  in 
determining  the  total  volume,  and  the  intrapulmonary  distribution 
of  the  inspired  gas  reaching  the  alveolar  gas  exchange  surfaces. 

It  is  known  that  the  intrapulmonary  distribution  of  inspired  gas 
is  altered  by  such  factors  as  1)  body  position;  2)  lung  volume; 

3)  chest  wall  geometry  and  compliance;  and  4)  increased  gravita- 
tional-inertial force  environments. 

Because  of  the  very  fragile  anatomic  structure  of  the  lungs  and 
the  large  intrathoracic  pressure  differentials  generated  by 
changes  in  the  force  environment  due  to  the  extreme  differences 
in  specific  gravity  of  the  gaseous  contents  of  the  alveoli  and 
the  parenchyma  and  blood  surrounding  these  alveoli,  the  lungs  are 


the  most  susceptible  organ  system  in  the  body  to  malfunction  and 
structural  injury  due  to  acceleration.  Furthermore,  it  is  quite 
possible  that  the  regional  distribution  of  mechanical  stress  within 
the  lung  is  a predisposing  factor  in  the  localization  of  pulmonary 
pathology.  Therefore,  a quantitative  description  of  the  dynamic 
spatial  distribution  of  lung  parenchymal  stresses  and  strains  and 
the  associated  changes  in  dynamic  shape  and  dimensions  (i.e.,  the 
dynamic  geometry)  of  the  intact  lungs  and  thorax  is  required  before 
further  progress  can  be  made  in  describing  the  three-dimensional 
mechanical  behavior  of  the  intact  lung.  However,  elucidation  of 
the  intrinsic  elastic  behavior  of  the  lung  and  its  relationship 
to  the  chest  wall  and  diaphragm  require  studies  of  the  intact  lung 
with  techniques  possessing  sufficient  temporal  and  spatial  resolu- 
tion to  permit  quantitative  measurements  of  the  dynamic  changes  in 
the  three-dimensional  parenchymal  stresses  and  strains,  regional 
lung  volumes,  and  thoracic  shape  and  dimensions  along  with  simul- 
taneous measurements  of  transpulmonary  and  intravascular  pressures. 

To  date,  most  of  the  knowledge  of  the  mechanical  properties  of 
the  lung  in  the  intact  thorax  has  been  derived  from  pressure-volume 
curves  (i.e.,  "compliance"  measurements)  (1 , 2)  and  the  regional 
distribution  of  an  inspired  radioactive  gas  (e.g. , xenon-133 ) (1, 3 , 4 ) . 
However,  these  methods  do  not  possess  sufficient  temporal  and 
spatial  resolution  to  allow  study  of  regional  lung  displacements, 
strains  and  volumes.  Hence,  dynamic  studies  of  the  intact  lung 
have,  heretofore,  not  been  possible. 

We  have  developed  in  our  laboratory  a method  (5,6)  for  determining 
the  spatial  and  temporal  distributions  of  lung  parenchymal  strains 
and  volumes  within  the  intact  thorax  by  recording  the  dynamic 
displacements  of  radiopaque  metallic  markers  by  means  of  a biplane 
videoroentgenographic  recording  system  which  permits  high  temporal 
(60-per-second)  and  spatial  (+1.0  mm)  resolution  measurements  of 
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the  "tagged"  lungs  during  various  respiratory  maneuvers,  altera- 
tions in  posture,  and  changes  in  gravitational-inertial  force 
environments.  This  "parenchymal  marker"  technique  has  been  used 
to  study  regional  parenchymal  strains  in  anesthetized,  as  well 
as  awake,  spontaneously  breathing  dogs  and  results  from  these 
studies  to  date  are  summarized  below. 

The  technique  involves  tagging  the  lung  parenchyma  with  radio- 
paque markers  (5).  The  metallic  markers,  which  are  1-mm  in 
diameter,  are  implanted  percutaneously  under  fluoroscopic  control 
throughout  the  parenchyma  of  the  right  lung  of  anesthetized  dogs. 

A three-dimensional  grid-like  array  of  up  to  30  markers  is  obtained 
without  causing  physiologically  significant  pneumo-  or  hemothorax. 

Postmortem  pathologic  studies  have  revealed  no  pleural  adhesions 
following  implantation  of  the  markers  and  pulmonary  parenchymal 
pathologic  changes  are  limited  to  a region  100-200  micra  surrounding 
each  marker  (5) . A minimal  recovery  period  of  2-3  weeks  is  pro- 
vided before  the  dogs  are  used  in  an  experiment. 

The  dynamic  displacements  of  the  parenchymal  markers  are  measured 
by  means  of  a biplane  videoroentgenographic  system  and  recording 
equipment  shown  in  Figure  2.  The  two  x-ray  sources  can  be  arranged 
in  either  an  orthogonal  or  in  a stereo  arrangement  as  shown.  The 
x-ray  sources  are  pulsed  on  nearly  simultaneously  and  transradiate 
the  thorax.  The  respective  fluoroscopic  images,  recorded  on  two 
9-inch  image  intensif iers,  are  scanned  by  television  (Plumbicon) 
cameras  and  the  biplane  videoroentgenograms  are  combined  into  a 
single  biplane  video  format  (Figure  3)  by  special  electronic 
switching  and  delay  circuitry  (7) . 
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Figure  2 Schematic  drawing  of  biplane  videoroentgenographic 

system  used  to  study  dynamic  regional  lung  mechanics. 

The  sequence  of  60/s  biplane  images  are  recorded  on  videotape  in 
absolute  temporal  synchrony  with  up  to  16  channels  of  analog 
information,  such  as  intrathoracic  pressure,  ECG  and  ventilation 
at  sampling  rates  as  high  as  1000  samples/s  per  analog  channel  by 
means  of  a video  amplitude  modulator-demodulator  system  (8). 


Since  the  image  intensifier  causes  pincushion  distortion  in  the 
recorded  video  images,  which  affects  the  accuracy  of  the  geometric 
measurements  of  the  distribution  of  the  metallic  markers  within 
the  image,  a method  has  been  developed  to  express  the  magnitude 
and  direction  of  this  distortion  over  the  entire  image  field 
covered  by  the  face  of  the  9-inch  image-intens if ier  screen.  A 
fourth  order  polynomial  with  15  terms  relating  selected  locations 
in  the  image  to  their  true,  x-y  coordinates  is  used  to  correct 
for  the  pincushion  distortion  in  the  recorded  x-ray  video  images. 
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This  predetermined  polynomial  can  correct  the  geometric  distor- 
tions (ranging  from  zero  at  the  center  of  the  field  to  about 
1.5  cm  at  the  periphery  of  the  image)  to  less  than  1 mm  error 
for  the  entire  video  field.  Image  magnification  due  to  diver- 
gence of  the  x-ray  beam  as  it  traverses  from  the  finite-size 
source  to  the  image-intensif ier  screen  is  also  corrected  (5) . 

To  facilitate  analysis,  the  biplane  images  and  multichannel  analog 
recordings  are  transferred  from  videotape  to  stop-action  video  disc 
(Ampex  DR-10A) , and  the  tracking  of  the  geometric  positions  of  the 
metallic  markers  is  performed  with  the  aid  of  an  operator-interactive 
computer  program.  A computer-generated  cursor  is  displayed  on  a 
video  scan  converter  (Tektronix  4501)  which  "mixes"  the  cursor  with 
the  stop-action  display  of  the  biplane  thoracic  image.  The  inves- 
tigator moves  this  cursor  until  it  coincides  with  one  of  the 
metallic  markers  in  the  video  image  and  then  depresses  a key  on 
the  keyboard  which  indicates  to  the  computer  (CDC  3500)  that  the 
two-dimensional  geometric  coordinates  of  the  cursor  and  thus  the 
metallic  marker  are  to  be  read  and  stored  in  an  array.  This 
procedure  is  repeated  until  the  geometric  coordinates  of  all 
identifiable  markers  on  both  projections  of  the  lungs  for  a single 
video  frame  (1/3 0th  s)  are  input  into  the  computer,  at  which  time 
the  computer  also  digitizes  and  stores  the  up  to  16  channels  of 
analog  data  that  had  been  recorded  simultaneously  with  the  video 
information.  The  coordinates  of  each  marker  are  corrected  for  pin- 
cushion and  magnification  distortions  and  then  the  three- 

dimensional  spatial  coordinates  of  each  marker  within  the  thorax 
are  determined.  An  example  of  a video  frame  used  for  such  a pro- 
cedure is  shown  in  Figure  3. 


Photograph  of  single  stop-action  frame  of  stereo 
biplane  videoroentgenogram  illustrating  the  stereo 
biplane  view  of  the  canine  right  lung  with  metallic 
parenchymal  markers  in  both  images.  Computer- 
generated cursor  dots  used  to  input  geometric  posi- 
tions of  markers  shown  in  right-hand  image. 


The  simultaneously  recorded  analog  information  is  contained  in  the 
"ribbon"  at  the  normally  unused  left  margin  of  the  video  image 
while  identification  information  (e.g.,  time  and  code),  added  during 
the  experiment,  is  displayed  at  the  bottom  of  the  image.  In  this 
particular  photograph  the  geometric  positions  of  the  metallic  markers 
at  end  expiration  are  indicated  by  the  bright  computer-generated 
(cursor)  dots  in  the  right  panel  of  this  biplane  thoracic 
image.  Following  determination  of  these  points,  the  video  disc 
was  stepped  forward  to  end  inspiration  to  show  the  displacement  of 


the  markers  from  these  points.  By  comparing  the  positions  of  the 
markers  at  end  expiration  as  indicated  by  the  white  (cursor)  dots 
and  at  end  inspiration  as  indicated  by  the  markers  themselves,  a 
qualitative  measure  of  the  direction  and  magnitude  of  parenchymal 
displacement  occurring  during  one  respiratory  cycle  is  obtained. 

The  quantitative  accuracy  of  this  technique  is  indicated  in 
Figure  4,  which  compares  the  actual  measured  distances  between 
markers  embedded  in  a Lucite  plate  which  was  rotated  in  the  biplane 
x-ray  system  and  the  computer-determined  distances  obtained  from 
the  videoroentgenographic  images.  This  graph  indicates  that  the 
computer  analysis,  which  removes  pincushion  and  magnification 
distortion  and  calibrates  the  corrected  coordinates  to  real  dimen- 
sions is  accurate  to  within  +1.0  mm.  In  addition,  calculations  of 
strains  between  the  embedded  markers  (i.e.,  strain  = 0)  during 
rotation  resulted  in  a mean  (N  = 294)  strain  of  0.004%  + 0.002  (SEM) 
for  markers  separated  by  1,  but  not  more  than  5 cm  after  correction 
for  magnification  and  pincushion  distortions.  In  our  studies,  only 
marker  pairs  separated  by  1 cm,  but  not  more  than  5 cm,  are  used 
in  calculating  dynamic  regional  parenchymal  strains. 
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Figure  4 Comparison  of  actual  and  computer-calculated  distances 
between  pairs  of  radiopaque  markers  fixed  at  known 
distances  within  a Lucite  plate. 
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magnitude  of  regional  parenchymal  displacements  during  a respiratory 
cycle.  Analysis  of  the  dynamic  displacements  of  these  metallic 
parenchymal  markers  can  be  used  to  determine  the  distribution  of 
regional  strains  within  the  lung  of  the  intact  thorax  during  various 
respiratory  maneuvers  and  under  different  gravitational-inertial 
force  environments. 


Figure  6 indicates  the  anatomic  positions  of  three  "marker-pairs" 

at  functional  residual  capacity.  These  pairs  of  markers  were 

selected  for  analysis  of  strain  in  three  regions  of  the  lung, 

i.e.,  apex,  midlung,  and  base.  The  dynamic  changes  in  distance 

in  three-dimensional  space  between  each  pair  of  markers  represents 

the  strain  (Al/Lq,  where  Al  is  the  measured  change  in  distance 

between  two  markers  at  any  point  in  time  and  L is  the  measured 

o 

distance  between  these  same  two  markers  at  end  expiration)  in  that 
region  of  the  lung  for  that  particular  orientation  of  the  marker  pair. 


(Dog,  14kg,  Morphine  - Pentobarbital  Anesthesia, 
Water-Immersion  Respirator, 

Head -Up,  12  breaths  /minute  ) 
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Figure  6 Location  and  orientation  of  3 pairs  of  markers  at 
FRC  selected  to  measure  strain  in  three  regions 
of  the  lung,  i.e.,  apex,  midlung,  and  base. 
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Distributions  of  regional  lung  parenchymal  strains  during  three 
respiratory  cycles  starting  at  approximately  FRC  (closed  symbols) 
and  during  three  respiratory  cycles  starting  at  FRC  + 200  ml  (open 
symbols)  are  shown  in  Figure  7. 


(Dog,  14kg,  Morphine  - Pentobarbital  Anesthesia,  Water  - Immersion 
Respirator,  Head-Up,  12  breaths  / minute  ) 
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Figure  7 Distributions  of  regional  lung  parenchymal  strains 
and  displacements  during  three  respiratory  cycles 
starting  at  FRC  (closed  symbols)  and  three  respira- 
tory cycles  starting  at  FRC  + 200  ml  (open  symbols) . 
Upper  three  panels  indicate  strain  distributions 
for  the  three  pairs  of  markers  shown  in  Figure  6 
in  the  apex,  midlung,  and  base,  respectively. 

Lower  three  panels  represent  absolute  distance 
between  marker  pairs  at  the  three  specified  lung 
levels . 


The  top  three  panels  indicate  strain  distributions  for  the  three 
pairs  of  beads  in  the  apex,  midlung,  and  base,  respectively.  It 
can  be  seen  that  parenchymal  strain  was  larger  at  the  lower  lung 
volume  than  at  the  higher  lung  volume.  The  lower  three  panels 
represent  the  absolute  distance  between  the  various  markers  during 


i 


-14- 


the  respiratory  cycle.  While  the  absolute  distance  between  marker 
pairs  was  increased  at  the  higher  lung  volume,  the  absolute  excur- 
sion or  change  in  intermarker  distance  during  respiration  was  less 
at  the  higher  lung  volume  than  at  the  lower  lung  volume.  The 
reproducibility  of  the  measurements  of  regional  distributions  of 
strain  from  cycle  to  cycle  can  also  be  seen  in  this  figure. 

The  anisotropy  of  parenchymal  strains  in  three  separate  lung  regions 
is  shown  in  Figure  8.  The  upper  two  panels  represent  the  antero- 
posterior and  the  lateral  projections  of  three  sets  of  three  markers 
each.  The  three  sets  of  markers  represent  regions  in  the  apex, 
the  midlung,  and  the  base  and  were  chosen  to  include,  for  each 
region,  pairs  of  markers  oriented  in  virtually  orthogonal  direc- 
tions. The  lower  three  panels  represent  the  computed  strain  for 
the  pairs  of  markers  shown  in  the  upper  panel.  These  three  panels 
represent  strains  for  two  marker  pairs  in  the  apex,  the  midlung, 
and  the  base,  respectively.  These  graphs  indicate  that  for  three 
sequential  breaths  the  marker  pairs  in  both  the  apex  and  the  midlung 
regions  exhibited  different  strains  relative  to  their  orientation 
within  the  lung.  In  the  base  of  the  lung  a paradoxical  movement  is 
exhibited  by  beads  15  and  18  where  a decrease  in  separation  occurred 
during  inspiration.  This  behavior  has  been  observed  in  other  parts 
of  the  lung  as  well  and  was  predicted  by  West  and  Matthews  (9)  for 
the  apical  region  but  not  for  the  basal  region  of  the  lung  of  a 
head-up  animal.  Thus,  this  technique  may  provide  the  experimental 
data  necessary  for  evaluation  of  the  theoretical  analyses  (9)  of 
stresses  and  strains  in  the  lung  caused  by  its  weight. 
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( Dog , 14  kg,  Morphine- Pentobarbital , Anesthesia,  Wa(er  - Immersion 
Respirator,  Head-Up,  12  breaths/minute) 


SECONDS 


Figure  8 Anisotropy  of  parenchymal  strains  in  three 

different  regions  of  the  lung.  Upper  two  panels 
show  location  of  three  sets  of  3 markers  chosen 
to  include  in  each  of  three  regions  (i.e.,  apex, 
midlung,  and  base)  pairs  of  markers  oriented  in 
virtually  orthogonal  directions.  Lower  three 
panels  represent  strains  for  pairs  of  markers 
in  upper  panel. 


Since  the  above  studies  were  performed  on  anesthetized  dogs  sup- 
ported head-up  in  a water-immersion  respirator,  the  question  as 
to  whether  the  gravitational  orientation  (i.e.,  head-up)  of  the 
dogs  as  well  as  the  presence  of  a column  of  water  surrounding  the 
animals  influenced  the  results  obtained  had  to  be  answered.  To 
do  this,  6 dogs  were  selected  for  training  prior  to  being  studied 
in  the  awake  standing  position.  The  parenchyma  of  the  right  lung 
of  these  dogs  was  tagged  with  radiopaque  metallic  markers  and  they 
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were  allowed  a minimum  of  3 weeks  to  recover.  After  this  period 
of  time,  the  animals  were  brought  to  the  laboratory,  on  an  indi- 
vidual basis,  and  allowed  to  run  free  for  a short  time  to  famil- 
iarize themselves  with  the  surroundings.  They  were  then  placed 
in  a Pavlovian-type  sling  (Figure  9)  and  positioned  on  the  fluoros- 
copic table  where  they  eventually  learned  to  stand  quietly  for 
20-30  minutes. 


Figure  9 Awake  dog  standing  in  a Pavlovian-type  sling  in 
front  of  the  biplane  videoroentgenographic  re- 
cording system.  Dog  can  also  be  anesthetized  and 
supported  in  this  "prone"  position  to  maintain  the 
same  orientation  of  the  lung  and  intrathoracic 
pressures  in  the  gravitational  field  as  in  the 
awake  state. 
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This  regimen  was  followed  for  3 consecutive  days,  but  on  the 
third  day,  the  x-ray  units  were  activated  and  biplane  videoroent- 
genographic  recordings  were  made  of  the  displacements  of  the 
parenchymal  markers  during  spontaneous  breathing. 

No  pressures  or  volumes  were  obtained  in  these  awake  animals  and 
no  attempt  was  made  to  influence  their  breathing  pattern  or  fre- 
quency. In  5 of  the  6 animals,  several  5-10  minute  segments  of 
spontaneous  breathing  were  recorded  in  each  animal  before  the 
animal  moved  and  had  to  be  repositioned. 

Continuous  tracking  of  the  displacement  of  individual  parenchymal 
markers  during  inspiration  and  expiration  revealed  that  for  the 
most  part  the  markers  returned  to  the  position  from  which  they 
started,  although  minor  degrees  of  irreversibility  were  seen  and 
these  were  corrected  within  2-3  breaths.  Occasionally,  a repro- 
ducible hysteresis  was  observed  in  the  displacement  pathways  of 
some  markers  during  inspirations  and  expiration.  In  addition, 
parenchymal  markers  located  near  the  heart  undergo  additional 
oscillation  and  displacement  due  to  cardiogenic  motion  with  this 
movement  being  particularly  visible  during  expiratory  pauses.  The 
extent  to  which  cardiac  motion  alters  the  stress-strain  relation- 
ship of  lung  parenchyma  adjacent  to  the  heart  has  not  yet  been 
determined.  Regional  lung  strains  during  spontaneous  breathing 
in  an  awake  standing  dog  are  shown  in  Figure  10. 


Markers  were  selected  to  measure  regional  parenchymal  strain  along 
the  horizontal  or  Z-axis  of  the  lung.  Markers  2,3,7,10  and  12  all 
lie  at  approximately  the  same  level,  that  is  on  the  same  isohydro- 
static pressure  plane  and  extend  over  the  major  portion  of  the 
apex-to-base  dimension  of  the  lung.  Note  that  the  strains  between 
successive  marker-pairs  along  this  line  are  very  similar  during 
the  two  breaths  plotted  here.  During  the  time  of  recording,  this 
animal  sighed  spontaneously  and  those  results  are  shown  in  Figure  11. 
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Figure  10 


DYNAMIC  REGIONAL  LUNG  STRAINS  IN  AWAKE  DOG 
(Dog.  13  kg.  No  Medication,  Standing  Position.  20-30  breaths  / minute  I 
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Dynamic  regional  parenchymal  strains  between 
markers  lying  along  horizontal  Z-axis  (i.e., 
isohydrostatic  pressure  plane)  of  the  canine 
lung. 
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Figure  11  Regional  lung  parenchymal  strains  between  markers 
lying  along  horizontal  Z-axis  of  lung  during 
spontaneous  sigh  in  awake  standing  dog.  Outlines 
of  lung  at  end-expiration,  taken  as  functional 
residual  capacity  (FRC)  for  this  animal  and  at 
total  lung  capacity  (TLC ) are  shown  to  indicate 
extent  of  lung  parenchyma  displacement  during 
this  particular  respiratory  maneuver.  Regional 
lung  parenchymal  strains,  determined  by  measuring 
the  dynamic  change  in  distance  in  three-dimensional 
space  between  marker-pairs,  is  shown  in  the  right- 
hand  panel. 


The  outlines  of  the  lung  at  end-expiration,  taken  as  FRC,  and  at 
the  peak  of  the  sigh  are  shown  to  indicate  the  extent  of  lung 
parenchyma  displacement  during  this  particular  maneuver.  Note 
again  the  near  equality  of  regional  parenchymal  strains  suggesting 
that  the  lung  behaves  like  a linear  spring  in  the  isohydrostatic 
pressure  plane  regardless  of  the  types  of  voluntary  respiratory 
maneuver  performed. 


The  anisotropy  of  regional  parenchymal  strains  can  be  seen  in 
Figure  12.  The  upper  two  panels  represent  the  dorsal-ventral 
and  lateral  projections  of  two  sets  of  3 markers  each;  namely, 

3,  6 & 7 and  7,  10  & 11.  These  two  sets  of  markers  were  selected 
to  represent  the  midlung  and  the  base  and  were  chosen  to  include 
in  these  2 regions,  pairs  of  markers  oriented  in  virtually  orthogonal 
directions.  The  middle  two  panels  represent  the  strains  between 
these  markers  during  3 consecutive  breaths.  Note  the  reproducibility 
of  the  regional  strains  from  cycle  to  cycle.  Note  also  that  strains 
between  marker-pairs  oriented  in  orthogonal  directions  (e.g.,  7,  11 
vs  7,  10  (right-hand  panel))  differ  by  75-80%,  with  these  differences 
being  significant  at  the  1%  level. 

Thus,  in  both  1)  anesthetized  dogs  supported  head-up  in  a water- 
immersion  respirator  and  in  2)  awake  dogs  in  the  standing  position, 
significant  anisotropy  of  regional  lung  parenchymal  strains  was 
observed.  In  other  words,  regional  lung  expansion  and  contraction 
is  not  necessarily  isotropic  or  uniform. 

This  anisotropic  behavior  reflects,  at  least  in  part,  the  influence 
of  the  shape  of  the  lung,  local  distortion  produced  by  the  beating 
heart,  diaphragm  and  chest  wall,  regional  variation  in  the  degree 
of  structural  inhomogeneity,  and  a non-uniform  spatial  distribution 
of  stress. 

Since  it  has  been  shown  that  the  mechanical  behavior  of  the  intact 
lung  is  dependent  on  absolute  lung  volume  and  the  lung's  previous 
volume  history,  it  is  necessary  to  measure  lung  volume  and  be  able 
to  relate  regional  mechanical  behavior  to  lung  volume.  To  do  this, 
we  have  designed  and  constructed  a radiolucent  lucite  air-filled 
plethysmograph  (Figure  13)  which  permits  quantitation  of  lung 
volume  and  can  be  used  with  anesthetized  spontaneously  breathing 
dogs  in  various  positions  (i.e.,  prone,  supine,  right  and  left 


lateral  decubitus)  and  with  awake  dogs  trained  to  lie  quietly 
with  the  plethysmograph  closed. 


(Dog.  8kg.  Stondmg  Position,  Promazine  2mg/kg.  i.m.) 


LUNG  OUTLINE 

at  FRC 

— at  Peak  of  Sigh 


Right  Left 

DORSAL  - VENTRAL 


£ CONDS 


Anisotropy  of  lung  parenchymal  strains  during 
spontaneous  breathing  and  a deep  sigh  in  an  awake 
standing  dog.  Upper  two  panels  show  locations  of 
markers  selected  for  analysis.  Lower  two  panels 
represent  strains  during  spontaneous  breathing 
(middle  row)  and  during  a spontaneous  sigh  (lower 
row)  . 


Figure  13  Radiolucent  plethysmograph  positioned  in  front 
of  the  biplane  videoroentgenographic  system. 

An  anesthetized  dog  is  shown  in  the  plethysmo- 
graph in  the  prone  position.  Plethysmograph 
permits  quantitation  of  lung  volume. 

Studies  with  anesthetized  dogs  in  different  positions  are  currently 
in  progress.  A measure  of  the  transpulmonary  pressure  (TPP)  or 
distending  pressure  of  the  lung  is  obtained  by  subtracting  airway 
pressure  from  pleural  pressure,  the  latter  being  provided  by  the 
pressure  recorded  from  an  esophageal  balloon  positioned  at  the 
mid-lung  level  (10,11). 
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The  animals  breathe  through  a pneumotachograph  in  order  to  measure 
velocity  of  flow  and  volume  (integral  of  flow)  of  air  moved  during 
various  respiratory  maneuvers.  Selected  breathing  regimens  (i.e., 
various  tidal  volumes  and  cycle  lengths,  including  periodic  respira- 
tory pauses)  and  volume  histories  are  imposed  and  controlled  by  the 
respirator  (Harvard)  during  positive  pressure  ventilation  (PPV) . 
Static  pressure-volume  curves  are  obtained  for  the  intact  lungs  by, 
starting  from  functional  residual  capacity  (FRC ) , inflating  the 
lung  with  a Hamilton  Super  Syringe  in  prescribed  increments  (e.g. , 
100  ml  increments)  to  total  lung  capacity  (TLC) , defined  as  25  cm 
H20  TPP  and  then  deflating  in  the  same  incremental  manner  back  to 
FRC.  Care  is  taken  to  impose  the  same  volume  history  on  the  lungs 
immediately  prior  to  doing  the  pressure-volume  curves. 


In  addition  to  using  pairs  of  markers  to  quantitate  directional 
inter-  and  intralobar  strains,  it  is  possible  to  calculate  dynamic 
changes  in  regional  lung  volume  (i.e.,  ventilation).  By  selecting 
sets  of  4 markers  (i.e.,  with  any  given  set  of  4 contained  within 
the  same  lobe)  to  define  a tetrahedron,  the  spatial  and  temporal 
distributions  of  changes  in  regional  lung  volume,  relative  to  the 
lung  itself  (i.e.,  in  contrast  to  inspired  radioactive  gas  tech- 
niques) are  obtained  throughout  the  respiratory  maneuver.  Further- 
more, since  the  three-dimensional  (i.e.,  spatial)  coordinates  of 
the  tetrahedrons  are  known  throughout  the  respiratory  cycle,  the 
dynamic  changes  in  relative  position  of  the  markers,  and  hence  the 
volume  described  by  them,  can  be  analyzed  with  established  methods 
in  engineering  mechanics.  For  example,  pure  deformation  can  be 
separated  from  rigid  body  motion  (i.e.,  translation  and  rotation) 


and  analyzed  for  axial  strains  (i.e.,  e 


e , e ) and  shear  strains 
y ' z 


(i.e.,  £Xy,  e , eXz)  t^iere^y  providing  strain  tensors  for  different 
regions  of  the  lung  during  a particular  respiratory  maneuver. 
Software  to  perform  these  calculations  has  been  written  and  imple- 
mented on  our  CDC  3500  computer. 


To  date,  3 dogs  have  been  studied  in  this  plethysmograph  and 
the  results  will  be  reported  later  upon  completion  of  this  phase 
of  the  study. 

To  aid  in  interpreting  regional  parenchymal  strain  data,  it  would 
be  helpful  to  be  able  to  visualize  in  three  dimensions,  the  dynamic 
displacements  of  the  parenchymal  markers  throughout  a selected 
respiratory  maneuver.  While  true  three-dimensional  display  devices 
are  still  in  the  developmental  stage,  it  is  possible  to  obtain  a 
pseudo  3-D  or  perspective  display  of  marker  displacements  by  looking 
at  a 2-D  image  of  the  respective  marker  positions  as  the  image  is 
rotated  about  a given  axis.  The  software  required  to  accomplish 
this  has  been  written  by  personnel  in  our  laboratory  and  the  pre- 
liminary results  are  shown  in  Figure  14.  The  picture  depicts  the 
anterior-posterior  (i.e. , ventral-dorsal)  view  of  a dog's  thorax 
(Z-axis  vertical)  with  the  dog's  right  side  on  the  left  of  this 
picture.  The  white  dots  depict  the  positions  of  12  markers  at 
functional  residual  capacity  (FRC)(i.e.,  at  end-expiration)  while 
the  red  dots  represent  the  positions  of  the  respective  markers  at 
total  lung  capacity  (TLC ) (i.e. , maximal  lung  volume) . The  pathways 
that  the  markers  took  in  going  from  FRC  to  TLC  and  back  are  shown 
by  the  lines  connecting  the  white  (FRC)  and  red  (TLC)  dots.  Only 
7 of  the  15  data  points  for  this  respiratory  maneuver  are  shown 
for  each  marker  for  the  sake  of  clarity.  Markers  in  different 
lung  lobes  are  identified  by  the  color  of  the  lines  connecting 
the  dots  (i.e.,  upper  lobe  - yellow;  middle  lobe  - blue;  and 
lower  lobe  - orange) . 
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Figure  14  The  ventral-dorsal  view  (Z-axis  vertical)  of 
a dog's  thorax  (dog's  right  side  on  left  side 
of  picture)  with  the  marker  positions  identified 
at  FRC  (white  dots)  and  at  TLC  (red  dots)  with 
the  pathways  taken  by  the  individual  markers  as 
the  lung  was  inflated  from  FRC  to  TLC  and  back 
to  FRC  indicated  by  the  lines  connecting  the 
data  points  (dots).  Lobar  location  of  markers 
identified  by  color  of  connecting  lines  (i.e., 
upper  lobe  - yellow;  middle  lobe  - blue;  and 
lower  lobe  - orange) . 


Figure  15  represents  the  same  data  as  shown  in  Figure  14  but  rotated 
90°  to  the  right  (right-lateral) (left  panel)  and  90°  to  the  left 
(left-lateral)  (right  panel)  to  further  illustrate  the  capabilities 
of  this  display  algorithm.  In  fact,  the  image  can  be  rotated  about 
any  axis  in  any  desired  angular  increment.  Furthermore,  it  is 
possible  to  visualize  the  dynamic  displacements  of  these  markers 
during  the  respiratory  maneuver  at  any  desired  angle.  It  is  hoped 
that  these  types  of  display  options  will  aid  in  assessing  regional 


lung  parenchymal  distortions  and  strains  as  well  as  elucidating 
the  role  that  lung  fissures  and  lobar  slippage  play  in  influencing 
regional  lung  mechanics. 


Figure  15  The  right-lateral  (left  panel)  and  left-lateral 
(right  panel)  views  of  the  marker  positions  in 
the  dog's  thorax  in  Figure  14.  See  Figure  14' s 
legend  for  description  of  color  codes  used. 


Summary  and  Significance  (Part  I) 

The  development  in  our  laboratory  of  the  parenchymal  marker  technique 
(Part  I)  as  well  as  the  ability  to  reconstruct  the  three-dimensional 
shapes  and  dimensions  of  the  heart,  lungs,  chest  wall  and  diaphragm, 
in  the  intact  thorax  (Part  II)  presents  a unique  opportunity  to 
study  dynamic  regional  lung  mechanics  and  lung-chest  wall-diaphragm 
interaction  under  various  conditions  of  lung  volumes  and  body  posi- 
tion and  thereby  provides  the  methodology  to  study  the  effects  on 
regional  lung  function  of  alterations  in  these  parameters  induced 
by  changes  in  the  gravitational-inertial  force  environment.  The 
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results  described  in  Part  I represent  significant  new  information 
regarding  regional  mechanical  properties  of  the  lung  and  clearly 
indicate  the  need  to  study  the  intrinsic  elastic  behavior  of  the 
lung  and  its  relationship  to  the  chest  wall  and  diaphragm  in  the 
intact  thorax  with  techniques  possessing  sufficient  temporal  and 
spatial  resolution  to  permit  quantitative  determination  of  the 
dynamic  changes  in  spatial  (three-dimensional)  lung  parenchymal 
strains,  regional  lung  volumes  and  lung-chest  wall  geometries. 

Mechanical  studies  of  the  intact  lung  have,  heretofore,  not  been 
possible.  These  data,  which  cannot  be  obtained  by  other  methods, 
provide  a quantitative  measure  of  the  mechanical  properties  of  the 
lung  under  normal  physiologic  conditions  (i.e.,  an  intact  thorax 
and  normal  intrathoracic  pressures)  at  1G  and  will  preface  similar 
planned  studies  under  conditions  of  increased  or  decreased  gravita- 
tional-inertial force  environments  encountered  in  aerospace  flight. 

Because  of  the  very  large  differences  in  specific  gravity  of  the 
air  in  the  alveoli,  and  the  lung  parenchyma  and  its  mobile  blood 
content  plus  the  anatomically  fragile  nature  of  the  lungs,  the 
respiratory  system  is  the  most  susceptible  of  all  bodily  organs  to 
malfunction  and  actual  structural  damage  during  exposure  to 
changes  in  the  direction  and  magnitude  of  the  gravitational- 
inertial  force  environment. 

Consequently,  increased  knowledge  concerning  the  mechanical  charac- 
teristics and  dynamic  three-dimensional  geometry  of  the  lung  is  a 
basic  requirement  for  understanding  the  factors  which  determine 
man's  reactions  and  tolerance  to  short  and  long  term  changes  in 
the  force  environment  such  as  encountered  in  conventional  aircraft 
and  in  manned  space  flight. 


I!  * 
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PART  I I : Dynamic  Spatial  Reconstruction  Tomography  for  Quantitative 

Studies  of  Structure  and  Function  of  the  Heart,  Lungs, 
and  Circulation 


A.  Overview  and  Current  Status  of  X-Ray  Techniques  for  Reconstruc- 
tion and  Display  of  Moving  Organ  Systems 

Two  major  activities  in  our  laboratory  during  the  past  four 
years  have  opened  the  door  to  development  of  procedures  for  true, 
dynamic  three-dimensional  imaging  of  the  heart,  lungs,  and  circula- 
tion; namely,  1)  development  of  the  capability  for  high  quality 
recording  and  rapid  computer  digitization  of  multiplanar  x-ray  video 
projection  images  (8,12)  and  comprehensive  display  of  quantitative 
image  and  functional  data  derived  therefrom,  and  2)  computer  imple- 
mentation of  a fan  beam  cross-sectional  reconstruction  algorithm, 
called  ART  III  (13,14),  which  is  similar  in  concept  to  that  employed 
in  the  first  EMI  brain  scanner  but  which  works  satisfactorily  with 
a relatively  low  number  of  two-dimensional  x-ray  fan  beam  projec- 
tions recorded  with  our  video  fluoroscopic  system. 

The  capability  provided  by  the  video  image  digitizing  system  for 
obtaining  high  spatial  resolution  (up  to  2,000  samples,  each  0.1  mm 
apart,  on  each  of  250  lines  per  video  field,  each  line  1 mm  apart) 
and  high  temporal  resolution  (60  video  fields  per  second)  measure- 
ments of  the  changes  in  x-ray  densities  and  the  distribution  of  these 
densities  throughout  organs  such  as  the  heart  and  lungs  is  required 
for  determination  and  display  of  the  true  dynamic  cross-sectional 
geometry  of  these  organs. 


In  addition  to  the  algebraic  reconstruction  technique,  ART  III  (13), 
we  have  implemented  on  our  CDC  3500  computer  a convolution  recon- 
struction algorithm  (15)  which  also  operates  on  x-ray  fan  beam  multi- 
planar  projection  data.  If  enough  views  can  be  collected  (60  or 
more) , the  convolution  algorithm  produces  more  accurate  cross- 
sectional  images  in  less  time  than  are  obtained  with  the  algebraic 
technique.  For  64  x 64  reconstructions  from  60  views,  approximately 
3 minutes  are  required  if  ART  III  is  used,  whereas  about  1 minute 
is  required  if  the  fan  beam  convolution  algorithm  is  used. 

Figure  16  diagramatically  illustrates  the  system  developed  in  our 
laboratory  for  generation  and  collection  of  two-dimensional  multi- 
planar  x-ray  projection  images  for  dynamic  spatial  reconstructions 
of  the  thorax  using  a computer-controlled  single-source  single- 
detector x-ray  fluoroscopic  system.  The  x-ray  imaging  chain,  shown 
at  the  bottom  of  the  figure,  is  a new  high  performance  video- 
fluoroscopic  system,  unique  in  its  design  and  construction,  which 
is  called  the  SSDSR  for  single  source  dynamic  spatial  reconstructor 
(see  Appendix  I) . 

This  assembly  has  the  advantage  that  the  entire  thorax  of  medium 
sized  dogs  can  be  imaged  on  the  large  12"  x 12"  flat  screen 
(GdSC^)  with  minimal  geometric  distortion.  This  cannot  be  done 
with  the  conventional  9"  diameter,  curved  screen  electrostatic 
image-intensif ier  tubes.  In  addition,  the  two-stage  magnetically- 
focused  image-intensif ier  and  image-isocon  video  camera  have  far 
superior  dynamic  range,  spatial  resolution,  and  lag  characteristics 
by  comparison  with  conventional  fluoroscopic  systems  (16,  Appendix  I). 
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D I AG  RAM  OF  COM  PUT  E R - CONT  ROLLED  SINGLE  - SOURCE 
SINGLE  - DETECTOR  X-RAY  VIDEO  SYSTEM  FOR  COLLECTION 
OF  MULTIPLANAR  ROENTGEN  VIDEO  PROJECTIONS  FOR 
DYNAMIC  SPATIAL  RECONSTRUCTIONS  OF  THE 
INTACT  THORAX  AND  ITS  CONTENTS 
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Figure  16  Diagram  of  computer-controlled  single-source 
dynamic  spatial  reconstruction  system  (SSDSR) 


The  object  to  be  reconstructed  is  positioned  in  the  x-ray  field 
of  the  SSDSR  and  rotated  under  computer-control  about  an  axis 
perpendicular  to  the  central  plane  of  x-ray  beams.  The  fluoro- 
scopic image  at  each  successive  angle  of  view  is  scanned  by  the 
video  camera,  producing  an  electronic  video  signal  representation 
of  the  image,  which  is  made  up  of  up  to  250  scan  lines  for  maximum 
(60/sec)  temporal  resolution  reconstructions  of  the  beating  heart 
or  up  to  500  horizontal  scan  lines  for  reconstructions  of  slower 
moving  structures  (17). 


Each  line  in  the  video  image  represents  the  projections  of  an 
approximately  0 , 5 mm  thick  section  of  the  object  scanned.  Measure- 
ments of  the  transmitted  x-ray  intensities  on  each  video  line  are 
made  by  the  video  digitizer  (12)  and  are  repeated  on  the  same 
lines  for  each  successive  view  of  the  object  as  it  is  rotated  in 
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equi-angular  increments  through  a total  arc  of  180°  to  360°.  The 
video  digitizer  can  digitize  up  to  1000  points  on  a single  line  in 
real  time,  or  the  entire  video  image  .in  about  10  seconds.  These 
x-ray  transmission  measurements  provide  #the  data  required  for 
determination  of  up  to  500  parallel  adjacent  cross  sections  of  the 
object,  or,  in  other  words,  its  complete  three-dimensional  recon- 
struction. 

Reconstruction  of  a moving  object,  such  as  a beating  heart  within 
the  thorax  of  a live  animal,  using  a single  x-ray  source  imaging 
system,  requires  maintenance  of  exact  temporal  relationships 
between  physiological  events  and  x-ray  projection  image  recording. 
This  is  accomplished  by  slaving  all  components  of  the  system,  in- 
cluding the  computer,  to  the  master  oscillator,  labeled  SYNC  in 
Figure  16,  which  controls  the  video  sweep  circuits  of  the  video 
camera.  The  rate  and  phase  of  the  cardiac  and  respiratory  cycles 
are  kept  constant  and  synchronized  with  the  rotation,  x-ray  pulses, 
video  recording  and  digitizing  process  by  computer-controlled 
ventilation  of  the  lungs  and  pacing  of  the  atria  and  ventricles, 
all  in  exact  phasic  relationships  with  each  other  and  with  the 
60-per-second  oscillator  pulses.  The  ECG  and  other  hemodynamic 
variables  are  simultaneously  recorded  with  the  video  images  on  a 
video  tape  or  video  disc  (8)  so  that  sets  of  images  recorded  at 
the  same  phase  of  successive  respiratory  and  cardiac  cycles  at 
each  incremental  angle  of  view  can  be  selected  for  each  60-per- 
second  point  in  time  throughout  the  respiratory  and  cardiac  cycles 
to  obtain  dynamic  spatial  reconstructions  of  the  full  anatomic 
extent  of  the  thorax  and  its  contents. 

Figure  17  schematically  illustrates  how  a single  transverse  cross 
section  of  an  object  like  the  heart  may  be  determined  using  an 
algebraic  reconstruction  technique. 
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Cross-Section  Reconstruction  from  Divergent 
Geometry  Transmission  Absorption  Profiles 
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Figure  17  Diagram  of  geometric  considerations  for  trans- 
axial  x-ray  reconstruction  of  a cross  section 
of  the  heart. 

The  transaxial  plane  traversed  by  the  x-ray  beam  which  contains  the 
desired  cross  section  of  the  object  to  be  reconstructed  is  mathe- 
matically divided  into  a number  of  small  squares  (typically  from 
64  x 64  up  to  128  x 128)  and  the  intensity  of  the  x-ray  beam  trans- 
mitted through  that  plane  is  measured  at  several  points  (typically 
from  200  to  500)  along  the  detector  for  each  rotational  position 
of  the  x-ray  source  and  detector  with  respect  to  the  obiect.  The 
geometric  considerations  are  the  same  if  the  source  and  detector 
are  rotated  about  the  object  or  if  the  object  is  rotated  with 
respect  to  a fixed  source  and  detector.  The  transmitted  x-ray  bear 
measurements  from  the  detector  are  related  to  measurements  of  the 
incident  x-ray  beam  (obtained  by  removing  the  object  from  the  system) 
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for  all  angles  of  view  by  a system  of  linear  equations  in  which  the 
unknowns  are  the  x-ray  beam  attenuation  coefficients  in  each  small 
square  of  the  grid.  The  algebraic  solution  of  this  system  of 
equations  provides  an  array  of  numbers  which  are  proportional  to 
these  x-ray  attenuation  coefficients  and  which  can  be  displayed  as 
a matrix  of  picture  elements  with  individual  brightness  levels 
proportional  to  the  attenuation  coefficients  in  each  small  square. 
Such  a display  results  in  an  image  of  the  cross-sectional  structure 
of  the  object  scanned.  The  detector  in  this  figure  can  be  con- 
sidered to  be  1 video  line  of  the  up  to  500  video  lines  in  our 
x-ray  fluoroscopic  system,  each  of  which  can  be  used  to  reconstruct 
cross  sections  over  the  entire  anatomic  extent  of  the  heart  or 
thorax . 

The  left-hand  panel  in  Figure  18  is  a photograph  of  the  x-ray  video 
image  of  the  thorax  of  an  intact  dog  recorded  with  the  SSDSR  system 
during  rotation  of  the  dog  through  360°.  The  white  area  on  either 
side  of  the  chest  in  each  image  represents  the  unattenuated  x-ray 
beam  impinging  on  the  fluorescent  screen  and  demonstrates  the 
superior  dynamic  range  of  this  system,  by  comparison  to  conventional 
fluoroscopic  equipment,  inasmuch  as  the  spine  can  be  seen  through 
the  heart  shadow,  and  other  intrathoracic  structures,  including  the 
radiopaque  2 mm  diameter  catheters,  simultaneously  with  segments 
of  the  image  produced  by  the  unattenuated  x-ray  beam.  No  x-ray 
contrast  material  has  been  injected  into  the  dog  to  produce  this 
image. 
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Figure  18  X-ray  video  projection  image  of  intact  dog  thorax 
recorded  with  SSDSR  (left)  and  reconstructed  cross 
section  of  the  thorax  (right)  at  level  indicated 
by  brightened  line. 


The  vertical  column  of  bars  at  the  left  of  each  image  are  the 
amplitude-modulated  analog  signals  of  up  to  16  hemodynamic  variable 
including  ECG,  pressures,  etc.,  which  are  simultaneously  recorded 
with  the  video  images,  (8).  The  brightened  line  indicates  the  levej 
of  the  reconstructed  cross  section  of  the  thorax  shown  in  the  right 
panel . 


This  reconstruction,  produced  by  a fan-beam  convolution  algorithm, 
contains  117  x 117  picture  elements  and  was  determined  from  120 
views  recorded  over  a range  of  360°.  Note  the  spinal  vertebra  and 
spinal  cord,  ribs,  pleural  surfaces  of  the  lungs,  the  esophagus, 
bifurcated  trachea,  epicardial  surfaces  of  the  heart,  and  the  two 
catheters  (2-mm  diameter)  in  the  heart.  Some  trapped  air  intro- 
duced via  the  intracardiac  catheters  can  be  seen  in  the  chambers 
of  the  heart  of  this  dead  dog.  The  spatial  and  contrast  resolution 
achieved  in  this  reconstruction  is  very  encouraging,  considering 
that  conventional  (4  ma)  fluoroscopic  x-ray  levels  and  video  re- 
cording techniques  were  used. 


The  top  panel  in  Figure  19  depicts  12  brightened  0.5  nun  thick 
video  lines,  each  16 -nun  apart,  covering  the  anatomic  extent  of  the 
thorax  of  an  intact  dog.  These  levels  were  selected  from  one  360°, 
circumferential  scan  for  reconstruction  of  12  cross  sections  of  the 
dog's  chest.  The  video  lines  just  above  and  below  each  brightened 
level  were  also  digitized  and  averaged  with  the  selected  lines. 

Since  the  midpoints  of  adjacent  lines  are  spaced  at  1.0  mm  intervals 
averaging  of  3 lines  results  in  a cross  section  which  encompasses 
an  approximately  3 nun  thick  section  of  the  scanned  object. 


Figure  19  X-ray  video  projection  image  of  dog  thorax  with 
superimposed  lines  (top)  at  12  anatomic  levels 
selected  for  cross-section  reconstruction,  and 
12  separate  3-mm  thick  transverse  sections  through 
thorax  (bottom)  reconstructed  at  selected  levels, 
extending  from  the  apex  to  the  base  of  the  lungs. 
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The  bottom  panel  in  Figure  19  shows  these  12  reconstructed  cross 
sections  of  the  thorax  from  the  apex  of  the  lungs,  at  the  upper 
left  of  the  panel,  extending  cephalocaudally  to  the  base  of  the 
lungs,  at  the  lower  right  of  the  panel,  in  lef t-to-right , top-to- 
bottom  order.  Each  cross  section  contains  64  x 64  picture  elements 
and  was  reconstructed  in  approximately  1 minute  from  60  equally 
spaced  views  over  360°  by  a fan-beam  convolution  algorithm.  The 
spinal  vertebra,  ribs  as  they  course  through  the  plane  of  the  cross 
section,  the  esophagus  and  bifurcated  airway,  the  epicardial  sur- 
face of  the  heart,  the  pleural  surfaces  of  the  lungs,  and  the  dia- 
phragm can  all  be  visualized  at  the  appropriate  levels.  Air  in 
the  transverse  segment  of  the  large  intestine  can  be  seen  in  the 
lower  ro w of  cross  sections.  No  x-ray  contrast  material  was  in- 
jected into  the  dog  to  obtain  these  reconstructions. 

A very  useful  advantage  of  full  spatial  anatomic  reconstructions 
of  the  thorax  such  as  those  shown  in  Figure  19  is  that  multiple 
parallel  sections  oriented  at  any  desired  angle  in  relation  to  the 
x,  y and  z axes  of  the  thorax  can  be  computed  from  the  reconstructed 
transverse  images,  and  displayed  for  viewing  of  the  thorax  from  any 
desired  aspect. 

The  top  panel  in  Figure  20  is  a photograph  of  the  A-P  x-ray  video 
projection  image  of  the  dog's  thorax,  which  contains  all  of  the 
structures  of  the  thorax  superimposed  upon  the  recording  plane. 

The  bottom  panel  in  Figure  20  is  a computer-generated  display  of 
12  separate,  parallel  3-mm  thick  coronal  sections  through  the 
thorax,  extending  from  the  sternum  at  the  upper  left  to  the  back- 
bone at  the  lower  right.  These  3 mm  thick  coronal  sections  of 
64  x 64  picture  elements  each  were  computed  from  64  reconstructed 
transverse  sections  of  the  chest,  of  which  12  are  shown  in 
Figure  19.  Note  the  ribs  in  the  chest  wall,  the  catheters  in  the 
heart  (first  and  second  images  in  row  2),  the  trachea  and  esophagus 
(third  and  fourth  images  in  row  2,  respectively)  and  the  detail  of 
the  vertebral  column  (in  row  3) . 


Figure  2 0 X-ray  video  projection  image  of  dog  thorax  (top) 
recorded  in  anterior-posterior  position,  and  12 
separate  3-mm  thick  corona  1 sections  of  the  thorax 
(bottom)  extending  from  the  sternum  to  the  back- 
bone, computed  from  the  same  scan  data  and  spatial 
array  x-ray  attenuation  values  encompassinq  the 
full  thorax,  used  for  the  reconstructed  sections 
shown  in  Figure  19. 
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Figure  21  is  a similar  display  with  a lateral  x-ray  video  projection 
image  of  the  dog's  thorax  shown  in  the  top  panel  and  12  separate, 
parallel  3-mm  thick  sagittal  sections  shown  in  the  bottom  panel. 

These  sections  extend  from  the  left  side  of  the  thorax  (at  the 
upper  left)  to  the  right  side  of  the  thorax  (at  the  bottom  right) 
and  were  computed  from  the  same  set  of  64  reconstructed  parallel 
transverse  cross  sections  as  were  used  to  compute  the  coronal 
sections  shown  in  Figure  19.  The  intrathoracic  structures  are  again 
clearly  visible.  Note  particularly  the  shapes  of  the  heart  and 
lungs,  the  catheter  traversing  the  cephalo-caudad  extent  of  the 
superior  vena  cava  (2nd  image)  and  in  the  right  atrium  (3rd  and 
4th  images,  row  2),  and  the  detail  of  the  spinal  vertebrae,  also 
in  the  second  row  of  images. 

To  further  illustrate  the  advantages  of  full  spatial  (3-D)  recon- 
structions of  the  thorax.  Figure  22  shows,  in  the  top  panel, 

12  sections  through  the  cephalo-caudad  mid-line  (Z-axis)  of  the 
thorax,  with  each  successive  section  from  lef t-to-right  and  top- 
to-bottom  representing  a slice  directly  through  this  central 
vertical  axis  but  rotated  by  24°  with  respect  to  the  previous 
section.  The  12  images  shown  in  the  bottom  panel  of  Figure  22  are 
parallel  oblique  sections  computed  at  an  angle  of  45°  to  each  of 
the  x,  y,  and  z axes  of  the  thorax.  This  slice  orientation  is 
nearly  parallel  to  the  longitudinal  base-to-apex  axis  of  the  heart. 
The  multioriented  "thin"  sections  in  figures  19-22,  all  of  which  were 
computed  from  one  synchronous  cylindrical  scan  of  the  thorax, 
demonstrate  the  capability  for  detailed  quantitative  analysis  of 
simultaneous  anatomical  structure  and  functional  relationships  of 
intrathoracic  or  other  organ  systems  of  the  body. 
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X-ray  video  projection  image  of  dog  thorax  (top) 
recorded  in  left-lateral  position,  and  12  separate 
3-mm  thick  sagittal  sections  of  the  thorax  (bottom) 
extending  from  the  left  chest  wall  to  the  right 
chest  wall,  computed  from  the  same  scan  data  and 
spatial  array  of  x-ray  attenuation  values  deter- 
mined for  Figure  19. 
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Twelve  longitudinal  thoracic  sections  (top) 
computed  at  24°  angular  increments  through 
cephalo-caudal  midline  (z-axis)  of  dou  thorax, 
and  .12  parallel  oblique  thoracic  sections 
(bottom)  computed  at  45°  angle  to  x , y,  z axes 
of  thorax  from  the  same  scan  data  and  resultin 
x-ray  attenuation  values  as  used  for  Figure  19 


However,  in  order  to  accurately  study  dynamic  structural-functional 
relationships  of  moving  organ  systems,  it  is  necessary  to  obtain 
high  temporal  resolution  (i.e.,  60  to  100  times  per  second)  three- 
dimensional  reconstructions  of  these  structures.  This  can  be 
accomplished  by  the  SSDSR  in  dogs  during  imposed  conditions  of 
"physiological  stationar ity " produced  by  precise  computer  control 
of  the  rate  and  magnitude  of  the  respiratory  and  cardiac  cycles 
in  exact  interphasic  temporal  relationships  with  the  multiplanar 
x-ray  video  image  scanning  and  recording  procedures. 
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X-ray  video  projection  images  of  living  dog's 
thorax  (top  row)  recorded  in  anterior-posterior 
position  during  3 phases  of  respiratory  cycle, 
and  reconstructed  cross  sections  of  thorax 
(bottom  three  rows)  for  these  3 phases  of  respira- 
tory cycle  at  3 different  anatomic  levels  of 
thorax,  indicated  by  brightened  horizontal  lines. 
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Figure  23  illustrates  a sequence  of  thoracic  cross-sectional  recon- 
structions of  a living  dog's  chest,  with  no  x-ray  contrast  material 
injected,  obtained  from  temporally  correlating  sets  of  digitized 
multiplanar  x-ray  images  of  the  thorax  recorded  during  computer- 
controlled  ventilation,  cardiac  pacing,  and  rotation  of  the  dog  in 
angular  increments  of  6°  on  successive  respiratory  cycles  for  a 
total  of  31  cycles  (i.e.,  31  views  over  180°),  The  top  row  of 
projection  images,  from  left-to-right , were  recorded  at  the  end- 
expiratory,  mid-inspiratory,  and  end-inspiratory  phases  of  the 
respiratory  cycle,  respectively.  The  brightened  horizontal  lines 
on  each  image  depict  three  levels  near  the  apex,  center  and  base  of 
the  thorax  selected  for  reconstruction  at  these  three  phases  of  the 
respiratory  cycle.  Each  row  of  cross  sections  (reconstructed  using 
ART  III,  due  to  the  limited  number  of  views)  represents  the  three 
selected  anatomic  levels  of  the  chest.  The  top,  center,  and  bottom 
rows  correspond  respectively  to  the  apical,  mid,  and  basal  levels  of 
the  thorax.  Each  column  of  cross  sections  represents  the  three 
different  phases  of  the  respiratory  cycle.  The  left,  center,  and 
right  columns  corresponding  respectively  to  the  end-expiratory, 
mid-inspiratory,  and  end-inspiratory  phases  of  the  respiratory 
cycle.  Hence,  this  montage  of  images  illustrates  both  dynamic 
(in  the  horizontal  direction)  and  spatial  (in  the  vertical  direction) 
reconstructions  made  from  a single  circumferential  scanning  pro- 
cedure. 

The  top  two  panels  in  Figure  24  are  x-ray  video  projection  recordings 
of  a beating  heart  in  an  intact  living  dog  at  end-diastole  (left) 
and  end-systole  (right) . The  brightened  horizontal  line  depicts 
one  of  approximately  80  to  100  levels  of  the  heart,  extending  from 
base-to-apex , which  could  be  selected  for  cross-sectional  recon- 
struction of  the  heart  during  the  cardiac  cycle.  During  rotation 
of  the  dog  in  6°  increments  on  every  other  beat  through  a total 
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range  of  180° , yielding  31  views  in  a total  scan  time  of  approxi- 
mately 25  seconds  (heart  rate  was  150  beats/minut e ) , respiration 
was  suspended  and  x-ray  contrast  material  was  infused  continuously 
into  the  left  ventricular  chamber  to  achieve  continuous  opacification 
of  the  chamber  during  the  rotational  scan  period. 


Figure  24  X-ray  video  projection  images  of  beating  heart 
in  intact  dog  (top)  recorded  at  end-diastole 
and  end-systole  during  infusion  of  x-ray  contrast 
media  into  left  ventricle,  and  cross  section  of 
intact  beating  heart  (bottom)  reconstructed  at 
level  of  brightened  line  for  12  points  in  time 
throughout  the  cardiac  cycle. 
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The  lower  three  rows  of  panels  in  Figure  24  show  a single  cross- 
sectional  level  near  the  mid-plane  between  the  base  and  apex  of 
the  heart  which  was  reconstructed  (using  ART  III)  at  12  successive 
points  throughout  the  cardiac  cycle,  with  each  successive  point  in 
time  being  33  msec  apart. 

The  left  ventricular  chamber  and  epicardial  surfaces  of  the  heart 
can  be  seen  in  each  successive  image  in  the  first  two  rows  of  cross 
sections  during  the  isovolumic  and  contractile  phases  of  the  sys- 
tolic period,  followed  by  the  subsequent  diastolic  filling  period 
shown  in  the  bottom  row.  The  upper  left  hand  cross  section  corres- 
ponds to  end-diastole,  and  is  coincident  with  the  computer  pacing 
pulse  (or  "R"  wave)  which  initiates  the  contraction. 


These  cross  sections  demonstrate  the  heretofore  unobtainable  cap- 
ability for  obtaining  the  dynamic  cross-scctional  shape  and  dimen- 
sions of  both  the  epi-  and  endocardial  surfaces  of  the  intact 
working  heart  in  living  animals  over  the  entire  anatomic  and  tem- 
poral extents  of  the  myocardial  walls  and  cardiac  cycle,  respec- 
tively. 


In  order  to  study  left  ventricular  dynamics  in  greater  detail, 
ideal  roentgenographic  conditions,  using  an  isolated,  working 
canine  left  ventricle,  were  arranged.  The  top  row  of  panels  in 
Figure  25  are  x-ray  video  images  of  an  isolated  left  ventricle 
shown  at  end-diastole  and  at  end-systole  before  and  after  ligating 
the  anterior  branch  of  the  circumflex  coronary  artery.  The  two 
brightened  lines  near  the  mid-level  and  apex  of  the  isolated 
ventricle  depict  the  anatomic  levels  selected  for  reconstruction 
of  the  cross  sectionr  f the  ventricle  shown  in  the  bottom  two 
rows  of  panels.  The  coronary  circulation,  supplied  by  arterial 
blood  from  a donor  dog,  was  kept  completely  separate  from  an  arti- 
ficial systemic  circulation,  which  contained  10  Renovist  in  solution 
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with  isotonic  Ringer's  dextran  to  provide  continuous,  beat-to- 
beat,  roentgen  opacification  of  the  chamber.  During  computer- 
controlled  rotation  of  the  isolated  ventricle,  both  before  and 
after  the  circumflex  coronary  artery  was  tied  off,  the  heart  rate 
was  controlled  by  electrical  pacing  pulses  which  were  synchronized 
with  the  video  scan  rate  to  ensure  a fixed  temporal  relationship 
between  the  60-per-second  x-ray  images  recorded  at  each  angle  of 
view  and  the  mechanical  phases  of  each  successive  cardiac  cycle. 
The  ventricle  was  rotated  during  alternate  cardiac  cycles  so  that 
one  complete  cardiac  cycle  was  recorded  at  each  angle  of  view. 
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Figure  25  X-ray  video  projection  images  of  beating  isolated 
canine  left  ventricle  (top  row)  recorded  during 
diastole  and  systole  before  and  after  ligation  of 
circumflex  coronary  artery,  and  reconstructed 
cross  sections  of  ventricle  (bottom  two  rows)  at 
2 levels  indicated  by  brightened  lines,  determined 
for  diastole  and  systole  before  and  after  ligation 
of  the  circumflex  coronary  artery. 
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A1 though,  not  discernible  in  the  projection  images,  the  reconstructed 
cross-sectional  images  in  the  bottom  two  rows  of  Figure  25  indicate 
the  functional  effect  of  the  diminished  supply  of  blood  to  the 
anterior  wall  of  the  left  ventricular  myocardium.  In  the  first  row 
of  cross  sections,  which  show  the  mid-level  of  the  left  ventricle 
at  end-diastole  before  and  after  ligation  Cleft  two  panels,  respective- 
ly) and  at  end-systole  before  and  after  ligation  (right  two  panels, 
respectively) , the  anterior  and  posterior  papillary  muscles  can  be 
clearly  seen.  The  cross-sectional  area  of  the  chamber  at  this 
mid-level  is  noticeably  greater  in  the  images  reconstructed  after 
the  artery  was  ligated,  indicating  the  weakened  contractile  cap- 
ability of  the  ventricular  myocardium.  At  the  apical  level  (bottom 
row  of  cross  sections) , the  end-systolic  cross-sectional  area  is 
significantly  larger  in  the  image  reconstructed  after  ligation  of 
the  coronary  artery.  By  determination  of  the  borders  of  the  endo- 
cardial and  epicardial  surfaces  in  these  cross  sections  (by  tracing 
with  a light  pen  or  using  automatic  edge  detection  techniques  (21)) 
accurate  quantitative  analysis  of  regional  left  ventricular  structure 
and  function  (e.g.,  dynamic  chamber  dimensions,  wall  thickness, 
rate  of  wall  thickening,  contractile  patterns,  etc. ) can  be  carried 
out  and  studied  under  a variety  of  known  and  controlled  conditions. 

Figure  26  is  a computer-generated  gray  level  plot  of  two  cross 
sections  of  an  isolated  beating  left  ventricle  reconstructed  near 
the  base  (left)  and  apex  (right)  of  the  ventricle  during  the  same 
instant  in  the  cardiac  cycle.  The  myocardium  is  depicted  by  points 
with  different  gray  levels  which  are  proportional  to  the  degree  of 
stress  at  each  point  location  in  a direction  perpendicular  to  the 
plane  of  the  reconstructed  cross  section.  These  values  are  deter- 
mined by  mathematically  dividing  the  cross-sectional  image  of  the 
myocardium  into  a linked  set  of  several  hundred  triangular  elements 
and,  by  using  finite  element  analysis  techniques  (19),  computing  the 
stresses  at  each  node  of  each  triangle  on  the  basis  of  the  cross- 
sectional  shape  of  the  myocardium  and  the  simultaneously  measured 
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transmural  cardiac  pressures.  This  type  of  display  in  which  measur  J 
or  calculated  physiological  variables  are  superposed  upon  the  recon- 
structed cross-sectional  anatomy  permits  simultaneous  analysis  of 
structure  and  function  of  the  object  being  studied. 


Figure  26  Computer-generated  stress  plots  for  two  cross 
sections  of  left  ventricular  myocardium  recon- 
structed at  basal  (left)  and  apical  (right) 
levels  at  same  point  in  diastolic  phase  of 
cardiac  cycle. 


Display  of  the  three-dimensional  image  information  obtained  by 
computerized  tomographic  techniques  is  an  important  factor  in 
establishing  the  usefulness  of  these  capabilities  in  basic  research 
and  clinical  practice.  Time  delay  video  and  cine  recordings  of 
the  temporal  and  spatial  sequences  of  the  computer-gonerated 
cross-sectional  reconstructions  illustrated  in  these  figures  have 
been  made.  These  sequences  of  cross  sections  reconstructed 
over  the  full  anatomical  extent  of  the  beating  heart  and  breathing 


lungs  throughout  successive  cardiac  and  respiratory  cycles,  respect 
ly,  in  isolated  preparations  and  in  intact  dogs  can  be  viewed  in 
variable  time  base  modes  ranging  from  stop-action  to  real-time. 

Such  dynamic  displays  provide  the  capability  for  studying  the  spati 
relationships  and  instant-to-instant  changes  in  geometry  of  these 
organs  and  permits  evaluation  of  their  regional  and  integrated 
overall  function  in  the  intact  thorax  (20) . 


Examples  of  how  three-dimensional  images  can  be  displayed  to 
permit  analysis  of  the  spatial  relationships  and  regional  detail 
of  image  structures  is  illustrated  in  Figures  27  and  28. 


Figure  27  X-ray  video  projection  image  of  isolated  intact 
canine  heart  (top)  with  coronary  arteries  filled 
with  x-ray  contrast  material,  and  4 cross  sections 
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! of  heart  (center)  reconstructed  at  anatomic 

levels  indicated  by  brightened  lines,  and  perspec- 
tive 3-D  display  of  coronary  artery  tree  (bottom) 
computed  from  x,y,z  coordinate  positions  of 
coronary  arteries  in  48  reconstructed  cross 
sections  of  heart. 

The  top  panel  in  Figure  27  is  an  x-ray  video  image  of  an  isolated 
intact  canine  heart  whose  chambers  and  coronary  arteries  were 
filled  with  cotton  and  barium  paste,  respectively,  prior  to  rota- 
tion and  digitization  of  four  video  lines  selected  for  cross- 
sectional  reconstruction. 

The  center  panel  shows  the  reconstructed  cross  sections  at  the  four 
levels  indicated  by  the  brightened  horizontal  lines  in  the  top 
panel.  The  characteristic  cross-sectional  shapes  of  the  epi-  and 
endocardial  surfaces  of  both  the  left  and  right  ventricles  at 
each  level  are  accurately  reconstructed.  In  addition,  the  cross 
sections  of  the  coronary  arteries  as  they  traverse  the  myocardial 
walls  can  be  seen  as  bright  spots  (more  x-ray  dense)  at  each 
cross-sectional  level.  These  arteries  are  approximately  2-mm  in 
diameter  near  the  base  of  the  heart  (upper-left  cross  section) 
and  0.5  mm  in  diameter  near  the  apex  of  the  heart  (lower-right 
cross  section) . 

The  bottom  panel  in  Figure  27  is  a three-dimensional  perspective 
line  display  of  the  coronary  tree  of  the  heart  determined  from 
the  coordinate  positions  of  the  coronary  arteries,  which  were  entered 
into  the  computer  using  an  operator-interactive  display  and  light 
pen  assembly,  in  each  of  48  parallel  cross  sections  of  the  heart 
reconstructed  from  its  base  to  apex.  By  performing  rotational 
matrix  operations  on  these  sets  of  x,  y,  z coordinates,  line  re- 
presentations of  the  coronary  tree  can  be  viewed  from  any  orientation 
in  stop-action  or  during  continuous  rotation  to  facilitate  studies 
of  the  geometric  relationships  of  the  coronary  artery  system. 
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Similar  three-dimensional  displays  of  the  complete  coronary  vessels, 
obtained  by  entering  into  the  computer  complete  sets  of  reconstructed 
sections  encompassing  the  full  extent  of  the  large  coronary  arteries 
permits  computation  of  regional  cross-sectional  lumenal  areas  of 
these  arteries  for  more  accurate  detection  and  quantitation  of 
regional  stenosis  or  occlusions. 

Application  of  three-dimensional  boundary  surface  detection  and 
three-dimensional  computer  graphics  techniques  to  three- 
dimensional  reconstruction  data  (21)  provides  the  capability  of 
three-dimensional  displays  of  solid  objects. 

The  upper  left  panel  in  Figure  28  is  a computer-generated 
three-dimensional  surface  display  of  an  isolated  canine  heart 
which  was  determined  by  automatic  three-dimensional  boundary  detec- 
tion of  30  reconstructed  cross-sectional  levels  of  the  heart. 

The  upper  right  panel  in  Figure  28  is  a similar  display  of 
the  same  heart  rotated  90°  and  tilted  forward  15°.  The  bottom 
panel  in  Figure  28  is  a three-dimensional  gray-level  surface 
display  of  the  same  heart  but  mathematically  opened  into  halves 
to  permit  viewing  of  the  structural  detail  of  the  endocardial 
surfaces . 

The  capability  provided  by  the  video  image  digitizing  system 
(12)  for  obtaining  high  spatial  (up  to  2,000  samples  on  each 
of  250  lines  per  video  field)  and  high  temporal  (60  video  fields 
per  second)  resolution  measurements  for  cross-sectional  recon- 
struction of  moving  organs  also  permits  the  study  of  the  appearance, 
mean  transit  and  clearance  of  x-ray  contrast  substances  through, 
for  example,  the  myocardium  with  a degree  of  anatomical  and 
temporal  precision  not  possible  by  current  radioisotope  techniques. 


Figure  28  Computer-generated  three-dimensional  gray- level 
displays  of  isolated  intact  dog  heart  determined 
from  30  reconstructed  cross-sectional  images 
extending  from  base  to  apex  of  heart. 

This  capability  is  indicated  in  Figure  29  which  illustrates 
the  quantitative  determination  by  dynamic  background  subtraction 
of  the  spatial  distribution  of  the  coronary  circulation  within 
the  myocardium  of  a patient  undergoing  selective  coronary 
arteriography.  The  computer-generated  difference  images  (4  right- 
hand  panels)  of  contrast  distribution  at  the  same  specific  phase  ol 
four  successive  heart  beats  were  obtained  by  subtraction,  sample  for 
sample,  of  the  digitized  arrays  of  the  actual  video  images  (4 
left-hand  panels) , recorded  during  and  following  the  injection  of 
contrast  material,  from  the  digitized  arrays  of  the  video  images 


irouotatsa  juct.  |>i  ioi  to  Um  injection.  Simultaneous  ECG  recordinq 
on  the  videotape  allows  precise  temporal  and  phasic  registration 
between  the  pre-injection  video  images  and  the  post-injection 
video  images  which  are  subtracted  to  produce  difference  pictures 
at  each  point  in  the  cardiac  cycle  throughout  successive  post- 
injection cycles. 


Hear  t 
Cycle: 

I 


ANGIOGRAM  SUBTRACTION 


Figure  2lJ  X-ray  video  images  of  human  heart  (left  column) 

recorded  at  same  phase  of  diastole  in  1 successive 
cardiac  cycles  during  selective  coronary  angio- 
graphy, and  computer-generated  difference  images 
of  the  heart  (right  column)  determined  by  dynamic 
background  subtraction  for  this  same  phase  of  these 
4 cardiac  cycles  before,  during,  and  following  injec 
tion  of  x-ray  contrast  material  into  the  left  anteri 
descending  coronary  artery. 
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The  display  of  the  difference  images  greatly  enhances  the  ability 
to  see  the  myocardial  perfusion  patterns  of  the  microvascular 
anatomy  and  demonstrates  the  potential  capability  for  quantitative 
evaluation  of  myocardial  ischemia  and  infarction. 

By  combining  dynamic  three-dimensional  reconstruction  techniques 
with  dynamic  subtraction  techniques,  the  instant- to-ins tan t 
regional  changes  in  shape  and  dimensions  and  transmural  perfusion 
of  the  myocardium  (and  other  organ  systems)  can  be  accurately  and 
completely  studied. 

B.  Physiologic  Considerations  in  the  Trade-Off  Between  Temporal 
and  Spatial  Resolution  for  Reconstruction  of  the  Dynamic 
Cardiac,  Pulmonary  and  Circulatory  Anatomy  and 
Simultaneous  Functions 


An  ideal  computerized  transaxial  scanner  should  be  capable  of 
producing  reconstructions  of  multiple  synchronous  parallel  sections 
in  any  desired  orientation  and  encompassing  the  complete  anatomic 
extent  of  the  structure  of  interest.  These  juxtaposed  cross  sections 
should  have  sufficient  three-dimensional  spatial  and  temporal 
resolution  for  visualization  and  quantitative  study  of  anatomic  static 
and  dynamic  structural  relationships  to  biologic  functions  of  all 
macroscopic  anatomic  structures  of  interest  including  the  heart, 
lungs,  vascular  anatomy  and  circulatory  functional  dynamics  in  all 
regions  of  the  body. 

Since  reconstruction  images  are  projection  images,  they  do  present 
ambiguous  information  due  to  superposition  of  structures  or  non- 
uniformities in  the  axial  direction,  just  as  in  conventional  pro- 
jection images,  in  proportion  to  the  thickness  of  the  imaged 
structure. 
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The  spatial  and  density  resolution  and  precision  of  current  com- 
mercial scanners  is  impressive  (22) , but  only  for  stationary  objects 
with  concomitant  uniform  structure  in  the  axial  direction  and 
uniform  cross-sectional  dimensions  throughout  and  parallel  to  the 
axial  extent  of  the  cross  section.  Since  biologic  structures 
usually  exhibit  three-dimensional  structural  and  dimensional  varia- 
bility the  precision  of  their  cross-sectional  reconstruction  depends 
on  the  degree  of  this  variability.  The  significance  of  this  limi- 
tation depends  upon  the  nature  of  the  study,  i.e.,  the  geometry  of 
the  organ  relative  to  the  scan  thickness,  motion  of  the  organ,  and 
the  type  of  information  which  it  is  hoped  to  derive  from  the  cross- 
sectional  image. 


The  present  study  was  undertaken  to  study  the  possible  magnitude 
of  the  effects  of  these  factors  on  the  ability  to  measure  the  true 
shape  and  dimensions  of  biologic  structures.  The  cylindrical 
scanning,  single  source,  dynamic  spatial  reconstruction  system 
(SSDSR)  (see  Appendix  I)  was  used  to  provide  the  data  for  most 
of  the  results  reported  here.  Constrained  ART  was  used  to  recon- 
struct the  cross-sectional  images  from  profile  scans  obtained  at 
thirty-five  angles  of  view  equally  spaced  around  a circumferential 
scanning  range  of  180°. 


Effect  of  Thickness  of  the  Section  Scanned 

Figure  30  shows  a stop-action  display  from  a video  disc  of  a mono- 
plane video  roentgenogram  of  a urethane  cast  of  the  myocardium  of 
the  left  ventricle  of  a dog  with  the  chamber  filled  with  a dilute 
Renovist  solution  (25  mgm  iodine/ml  solution)  in  saline.  The 
horizontal  plane  containing  the  central  fan  beam  of  the  roentgen 
fluoroscopic  system  is  designated  by  the  brightened  horizontal 
line  midway  between  the  base  and  apex  of  the  cast.  The  urethane 
cast,  which  has  the  complex  shape  and  dimens. ions  and  a roentgen 


density  similar  to  the  left  ventricular  myocardium,  has  the  advan- 
tage that  it  can  be  physically  sliced  into  multiple  parallel  cross 
sections  at  the  same  levels  (beads  glued  to  surface  of  cast  identi 
these  levels)  selected  for  cross  section  reconstructions.  The 
slices  of  the  cast  then  serve  as  true  standards  against  which 
the  shape  and  dimensions  of  the  cast  determined  from  reconstructed 
cross  sections  can  be  compared. 

(Chamber  Filled  with  Renovisl  Solution,  25  mg  I/ml) 


iqure  30  A single  videoroentgenographic  view  of  urethane 

cast  of  canine  left  ventricle  with  chamber  filled 

with  dilute  solution  (25  mgm  iodine/ml)  of  Renovisu 
in  saline.  The  cast  war-  transf  i s with  l r 
diameter  aluminum  rod  for  mechanical  circumferenti *1 
mu i i ; planar  "cy 1 i ndr i ca ) " scann ing  foj  - • nsti  ic- 

tion  of  multiple  juxtaposed  cross  section:-.  Tin- 
dark  circular  objects  arc  1 mm  lead  beads,  at  taolu  i 
to  the  "epicardial"  surface  for  fiducial  points. 


Figure  31  illustrates  the  accuracy  of  reconstructed  cross  sections 

of  this  urethane  cast  of  the  left  ventricle,  as  determined  by 

tracing  with  an  operator-interactive  graph  pen  the  reconstructed 

cross-sectional  outlines  of  the  epi-  and  endocardial  surfaces  of 

the  left  ventricular  cast  at  21  levels  from  base-to-apex,  and 

comparing  the  cross-sectional  areas  computed  from  these  outlines 

with  those  computed  from  planimetry  measurements  of  the  epi-  and 

endocardial  surfaces  of  the  corresponding  actual  cross  sections  of 

the  cast.  The  mean  difference  between  outlines  (shapes)  of  the 

actual  and  reconstructed  cross  sections  for  all  levels  was  0.6  mm 

+_  .39  mm  for  the  epicardial  surfaces  and  0.83  mm  + 0.59  mm  for  the 

endocardial  surfaces.  An  actual  slice  of  the  cast  taken  from 

its  central  portion  and  a reconstructed  slice  determined  near  the 

same  level  (indicated  by  brightened  line  in  Figure  30)  are  shown  at 

the  top  of  Figure  31.  Correlation  of  the  areas  determined  from 

the  actual  and  reconstructed  cross  sections  were  very  good  as  shown 

in  the  regression  analysis  plot  at  the  bottom  of  Figure  31.  Accuracy 

of  chamber  volume  and  myocardial  mass  were  determined  to  be  within 

5%  of  the  true  values,  the  major  limitation  responsible  for  the 

2 

uncertainty  being  the  1 mm  size  of  the  individual  pixels  in  the 
array  which  constitutes  the  image  (23).  Similar  accuracy  of  recon- 
struction has  been  shown  using  other  test  objects  with  both  simple 
and  complex  geometric  configurations  (24). 

The  effect  of  the  thickness  of  the  scanning  plane  is  illustrated  in 
Figure  32  which  shows  comparisons  of  0.3  mm  thick  slices  (1  video 
scan  line)  and  10.2  mm  thick  slices  (17  video  scan  lines  averaged 
together)  at  three  levels;  10  mm  above  the  central  plane  (top  pair), 
at  the  central  plane  (center  pair) , and  20  mm  below  the  central 
plane  (bottom  pair)  of  the  ventricular  cast.  The  anatomic  sites 
of  the  different  levels  are  indicated  by  brightened  lines  on  the 
vidcoroentgenographic  projection  image  of  the  urethane  cast  on  the 
left.  At  the  upper  level  both  epicardial  and  endocardial  surfaces 


are  coincidentally  uniform  and  nearly  perpendicular  to  the  scan 
plane.  For  this  reason  accuracy  of  the  0.3  mm  thick  and  the 
10.2  mm  thick  slices  are  very  similar  at  this  level. 
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CROSS  SECTION 

Cast  Reconstructed 


Figure  31  Comparison  of  actual  and  reconstructed  cross- 
sectional  images  of  urethane  cast  of  left 
ventricle  (top)  and  of  epi-  and  endocardial 
areas  of  actual  and  reconstructed  cross  sections 
of  cast  determined  for  21  levels  from  base-to- 
apex  (bottom) . 
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EFFECTS  OF  AXIAL  VARIATIONS  OF  LEFT 
VENTRICULAR  SHAPE  AND  SCAN  SECTION  THICKNESS 
ON  CROSS-SECTIONAL  RECONSTRUCTIONS 


Figure  32  X-ray  video  projection  image  of  urethane  cast 
of  left  ventricle  with  three  superimposed 
brightened  lines  at  levels  selected  for  cross- 
sectional  reconstruction  (left)  and  two  sets 
of  reconstructed  cross-sectional  images 
determined  at  these  three  levels  for  scan 
thickness  of  0.3  mm  and  10.2  mm,  respectively 
(right) . 


However,  at  the  central  level  the  epicardial  border  becomes  less 
distinct  in  the  thicker  10.2  mm  slice.  This  effect  is  most  notice- 
able in  the  lower  right  and  upper  central  segments  of  the  endo- 
cardial surface.  Also  the  difference  between  the  position  of  the 
endocardial  border  in  the  lower  right  quadrant  of  the  0.3  mm  and 
10.2  nun  thick  slices  is  most  marked,  with  the  difference  approxi- 
mating one-half  the  wall  thickness.  In  the  region  20  mm  below  the 
central  plane  both  the  epicardial  and  endocardial  surfaces  are  not 
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perpendicular  to  the  scan  plane.  Therefore,  the  axial  projection 
of  a thick,  10.2  nun  and  a thin,  0.3  mm  slice  differ  significantly 
at  this  level.  The  endocardial  border  in  the  10.2  mm  thick  section 
is  considerably  "out  of  focus,"  making  determination  of  the  true 
shape  and  dimensions  of  the  left  ventricular  model  at  this  level 
impossible  when  using  such  a thick  scan  plane. 

These  results  are  particularly  important  for  cross-sectional  imaging 
of  the  heart  in  the  intact  chest  since  the  base-to-apex  axis  of  the 
ventricles  is  normally  at  an  oblique  angle  to  the  long  axis  of 
the  body  which  would  conventionally  be  perpendicular  to  the  scan- 
ning plane  (s)  of  a computerized  tomographic  scanner.  Plence,  the 
epicardial  and  endocardial  surfaces  of  the  heart  will  usually  not 
be  perpendicular  to  the  plane  of  reconstruction,  further  requiring 
that  transaxial  reconstruction  slices  be  as  thin  as  possible  (1  mm) 
to  permit  accurate  measurements  of  the  shape  and  dimensions  of  the 
heart.  In  addition,  simultaneous  multiple  parallel  cross  sections 
are  needed  to  determine  the  true  three-dimensional  shape  and  dimen- 
sions of  a moving  organ  such  as  the  heart. 

To  study  the  effect  in  reconstructed  images  due  to  motion  of  the 
object  during  the  circumferential  scanning  procedure,  a mathematical 
test  object  was  contrived.  The  simulated  test  object  was  a cylinder 
7.62  cm  in  diameter  with  groups  of  one,  two,  four,  and  eight  milli- 
meter triangular  serrations  (to  simulate  trabeculae  carnae  and 
papillary  muscles)  inserted  at  90°  intervals  around  the  perimeter. 
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Parallel  projection  data  was  mathematically  generated  for  the 
simulated  test  object  at  equiangular  increments  of  1.8°  over  a 
range  of  180°.  Random  quantum  noise  to  simulate  the  photon  sta- 
tistics of  the  SSDSR  was  introduced  into  the  projection  data  assuming 
an  incident  x-ray  count  of  10^  x-ray  photons  per  sampled  ray  (16) . 
These  mathematically  derived  projection  data  were  used  to  reconstruct 
the  test  object  shown  in  panel  A of  Figure  33  using  a parallel 
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convolution  algorithm  with  a 150  x 150  picture  element  matrix 
(pixel  size  of  0.7  mm  x 0.7  mm) . The  8 mm  (upper  right  quadrant 
of  Panel  A) , 4 mm  (upper  left  quadrant) , and  2 mm  (lower  left 
quadrant)  serrations  are  all  clearly  visualized  in  this  reconstruc- 
tion. The  location  of  the  1 mm  serrations  (lower  right  quadrant) 
can  just  be  seen  but  the  serrations  cannot  be  separated,  probably 
because  the  pixel  size  of  the  reconstruction  is  nearly  the  same  as 
the  size  of  the  serrations. 


(Mathematical  Simulations  ) 

NO  MOTION  ROTATION 

(7.62  cm  Diameter  ) ( 5 D - J . 6°  ) 


TRANSLATION 

( S D - / . 6 mm  ) 


RADIUS  CHANGE 

( SO  - 2.8  mm  ) 


Figure  33  Reconstructed  images  of  simulated  cylindrical 
test  object  showing  artifacts  caused  by  motion 
during  the  scanning  procedure.  Mathematically 
derived  projection  data  of  the  simulated 
7.62  cm  diameter  plexiglass  cylinder  with 
8,  2,  4,  and  1 mm  triangular  serrations  in 
equispaced  groups  around  the  perimeter  were 
used  to  obtain  reconstructions  without  motion 
(panel  A),  with  rotational  motion  (panel  B), 
with  translation  motion  (panel  C) , and  radial 
motion  (panel  D) . 
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This  image  (.panel  A)  reconstructed  from  mathematically  accurate 
projection  data  serves  as  a reference  against  which  the  other 
images  in  the  montage  can  Be  compared  for  image  degradation  due 
to  motion.  The  effect  of  a random  rotation  around  the  center  of 
the  cylinder  (as  occurs  in  the  heart  during  certain  phases  of  the 
cardiac  cycle)  during  the  scanning  procedure  is  illustrated  in 
panel  B.  This  effect  was  achieved  by  substituting,  at  each  angle 
of  view,  mathematically  derived  projection  data  randomly  selected 
to  be  at  a greater,  identical,  or  a lesser  angle  of  view  than  the 
correct  angle.  The  random  variations  in  the  angle  of  view  at  each 
viewing  position  were  normally  distributed  with  a mean  of  zero  and 
a standard  deviation  of  3.6°  resulting  in  a movement  of  up  to  2 mm 
at  the  periphery  of  the  cylinder.  The  smooth  round  perimeter  of  the 
reconstructed  test  object  is  very  well  defined  as  would  be  expected 
with  motion  tangential  to  this  surface.  However,  the  2 mm  and  4 mm 
serrations  are  obliterated  and  a faint  shadow  of  the  8 mm  serration 
is  just  visible.  With  no  prior  knowledge  of  the  shape  of  the  struc- 
ture in  this  region,  a true  determination  of  the  dimensions  of  the 
cylinder  would  be  difficult. 

The  effect  of  random  two-dimensional  translational  movements  (as 
occur  during  asynchronous  respiratory  and  cardiac  cycles)  of  the 
cylinder  in  the  plane  of  the  reconstruction  is  illustrated  in  panel 
C of  Figure  33.  This  effect  was  achieved  by  randomly  shifting 
the  projection  data  at  each  angle  of  view.  The  random  shift  was 
normally  distributed  with  a mean  of  zero  and  a standard  deviation 
of  1.6  mm.  This  caused  a significant  loss  of  definition  of  the 
smooth  surface  of  the  cylinder  and  failure  to  reproduce  the  2 mm 
and  4 mm  serrations. 

The  effects  of  random  radius  changes  (to  simulate  slight  changes  in 
ventricular  volume  during  breath  holding  or  due  to  slight  sinus 
arrhythmia) , is  shown  in  panel  D.  This  effect  was  achieved  by 
compressing  or  expanding  the  projection  data  with  scaling  about  the 
center  point.  The  magnitude  of  the  radius  change  was  normally 
distributed  with  a mean  of  zero  and  a standard  deviation  of  2.8  mm. 
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There  is  a significant  loss  of  definition  of  the  smooth  boundary  in 
the  reconstructed  image  although  the  8 mm  serrations  are  clearly 
visible.  The  4 mm  serrations  can  be  detected  but  are  not  well 
defined. 

These  studies  of  thickness  and  motion  artifacts  in  reconstructed 
cross-sectional  images  using  actual  and  simulated  test  objects 
indicate  that  the  single  or  dual  slice  scan  mode  of  all  curi'ently 
available  commercial  scanners  coupled  with  their  poor  temporal 
resolution  severely  limit  the  application  of  these  scanners  to  the 
study  of  the  structure  and  function  of  moving  organ  systems  such  as 
the  heart,  lungs  and  circulation.  Even  if  adequate  single  cross- 
sectional  images  can  be  generated  using  gating  techniques,  the 
difficulties  associated  with  accurate  positioning  of  the  scanner  at 
multiple  parallel  levels  in  relation  to  the  organs  under  study, 
especially  the  heart,  will  be  difficult  to  overcome  because  of 
motion  of  the  subject  or  motion  of  the  internal  organs  of  the  body. 
For  these  reasons,  high  temporal  and  axial  resolution  synchronous 
cylindrical  scanning  is  practically  a mandatory  requirement  for 
accurate  studies  of  the  structure  and  function  of  the  dynamic 
organs  of  the  body  (22,25). 

Summary  and  Significance  (Part  II) 


The  development  in  our  laboratory  of  a single  source  computerized 
x-ray  fluoroscopic  system  (i.c.,  SSDSR)  for  three-dimensional 
reconstruction  and  display  of  moving  organs,  particularly  the 
heart,  lungs  and  circulation,  provides  the  potential  for  studies 
in  intact  animals  and  humans  of  the  dynamic  anatomy  and  associated 
functions  of  moving  organs  in  ways  never  before  achievable.  This 
is  because  up  until  the  development  of  computerized  transaxial 
tomography  there  was  no  technique  capable  of  accurate  determination 
of  the  true  shape  and  dimensions  of  these  organs  in  intact  animals 
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and  humans.  Three-dimensional  reconstruction  techniques  can  be 
applied  to  intact  animals  and  humans  in  a variety  of  postures 
(e.g.,  prone,  supine,  right  and  left  lateral  decubitus  positions) 
and  under  conditions  of  increased  gravitational-inertial  stress 
for  heretofore  impossible  investigations  of  the  relationships  of 
the  dynamic  three-dimensional  structural  geometry  to  the  physical 
and  biochemical  functions  of  both  moving  and  stationary  organ 
systems. 
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The  three-dimensional  anatomic  distribution  of  tissues  of  differing 
densities  within  the  human  body  can  be  visualized  by  generating 
images  of  multiple  parallel  cross  sections  using  the  technique  of 
computer-assisted  synchronous  cylindrical  scanning  tomography 
(Johnson,  et  al. , 1974;  Robb,  et  al.,  1974).  Existing  computerized 
transaxial  tomographic  systems  scan  one  or  two  cross  sections  at  a 
time  and  have  other  operational  restrictions  which  severely  limit 
use  for  three-dimensional  reconstruction  of  moving  organs  or  structures, 
especially  the  heart,  lungs  and  the  circulation  in  all  regions  of 
the  body.  Primarily,  these  limitations  (Robinson  1975)  consist  of 
the  duration  of  scan  required  (generally  greater  than  5 seconds) 
per  cross  section  imaged,  thickness  of  the  cross  section  (often 
greater  than  10  mm)  and  need  to  mechanically  move  either  the  patient 
or  the  scanning  assembly  to  image  parallel  cross  sections.  These 
constraints  preclude  effective  "stop-action"  (<0.01  sec)  imaging  of 
a three-dimensional  volume  containing  a moving  organ  such  as  the 
heart,  lungs,  or  flowing  blood.  Optimal  reconstruction  of  the 
intact  thorax  requires  imaging  of  its  full  anatomic  axial  and  trans- 
verse extents.  Consequently,  the  s^ts  of  x-ray  projection  "pro- 
files" at  each  angle  of  view  should  cover  an  area  of  at  least  30  x 
30  cm  to  include  all  or  a major  portion  of  most  chests.  In  addition 
and  ideally,  the  roentgen  density  values  of  the  projection  images 
recorded  at  each  angle  of  view  should  have  adequate  sensitivity  and 
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signal-to-noise  characteristics  to  permit  discrimination  of  roentgen 
intensity  differences  caused  by  a change  of  x-ray  beam  attenuation 
of  0.1%  or  less  (see  Ritmnn,  et  al.,  this  volume). 

A prototype  x-ray  imaging  system  has  been  designed  and  fabricated 
to  achieve  the  imaging  specifications  required  for  reconstruction 
of  the  three-dimensional  distribution  of  roentgen  densities  with- 
in all  or  a major  portion  of  the  transverse  and  cephalocaudal 
extents  of  the  thorax  and  at  frequent  intervals  of  time  throughout 
individual  respiratory  and/or  cardiac  cycles.  The  system  generates 
convcntiona 1 roentgen  video  fluoroscopic  images  for  direct 
visualization  and  recording  as  well  as  direct  high  fidelity  on-line 
analog-to-digital  conversion  of  a raster  of  up  to  1000  discrete 
values  covering  the  horizontal  extent  of  each  of  a maximum  of  250 
lines  per  video  field  (Gilbert,  et  al.,  1976;  Robb,  et  al. , 1973) 
at  each  "angle  of  view"  around  the  subject.  From  these  data  up  to 
250  parallel  adjacent  cross-sectional  images  can  be  computed.  The 
thickness  and  spacing  interval  of  each  cross  section  is  about  0.4 
and  0.9  run,  respectively,  with  a pixel  resolution  of  250  x 250  at 
each  of  the  250  levels  (Johnson,  et  al.,  1974). 

Each  of  these  62,500  pixels  defines  a finite  area  element  on  the 
surface  of  each  cross  section.  Since  each  cross  section  is  of 
finite  thickness,  each  pixel  also  represents  a finite  volume  elci  , 
i.e.,  a voxel,  in  each  cross  section.  The  total  three-dimensional 
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reconstructed  array  of  voxel  x-ray  density  values  of,  for  example 
the  entire  thorax,  provides  the  capability,  without  sequential  scan- 
ning procedures,  of  computing  sectional  images  with  as  many  dif- 
ferent orientations  as  desired  in  relation  to  the  body  axes  (e.g. , 
sagittal,  coronal  or  other  sectional  orientations)  in  addition 
to  the  transaxial  sections  reconstructed  from  the  original  scan 
(Robb,  et  al.,  1976)).  The  needs  and  requirements  of  an  all  elect- 
ronic cylindrical  scanning  assembly  to  acquire  all  the  necessary 
data  within  0.01  sec  and  with  a maximum  repetition  rate  of  60  cy- 
lindrical scans  per  second  have  been  described  elsewhere  (Ritman, 
et  al . , 1975 ) . 


METHOD 


A schematic  of  the  prototype  imaging  system  is  provided  in  Fig.  1. 
The  x-ray  pulses  generated  by  a conventional  rotating-anode  x-ray 
tube  either  pass  through  the  subject  and  are  attenuated  according 
to  the  integrated  roentgen  density  through  the  subject  or  pass 
alongside  the  subject  remaining  essentially  unattenuated.  The 
resulting  x-ray  image  is  projected  onto  a 30  cm  x 30  cm  flat 
fluorescent  screen  after  passing  through  a focused,  85  lines  per  mm 
15:1,  collimating  grid.  The  fluorescent  screen  is  of  the  gado- 


linium dioxy sulphide  terbium  activated  (Gd202S:Tb)  type  selected 
for  low  lag,  low  noise  and  high  x-ray- to-light  conversion  efficiency. 
In  addition,  the  light  output  is  predominantly  monochromatic  at 
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approximately  540  nanometers  which  is  particularly  convenient  in 
objective  lens  design. 


An  f/0.95  50  mm  focal  length  objective  lens  (f/0.75  50  mm  lens 
will  be  utilized  when  it  becomes  available)  is  used  to  focus  the 
fluorescent  screen  image  onto  the  input  photoelectric  surface  of 
a magnetically  focused  image  intensified.  The  objective  lens 
minifies  the  image  to  produce  a 28  x 28  mm  image.  In  an  assembly 
in  which  the  objective  lens  and  image  intensifier  are  positioned  at 
right  angles  to  the  fluorescent  screen,  a dichroic  mirror  was  found 
to  be  the  only  suitable  mirror.  Optical  coating  of  all  objective  lens 
surfaces  was  required  to  reduce  scatter  and  maximize  transmission.  In 
addition  the  inner  surface  of  the  structural  enclosure  housing  the 
objective  lens  and  mirror  had  to  be  baffled  with  light  traps  and 
painted  with  highly  absorbent  paint  (Nextel  Velvet  Coating,  3M)  to 
reduce  the  noise  produced  by  light  reflected  back  from  the  enclosure 
onto  the  fluorescent  screen. 


The  image  intensifier  amplifies  light  intensity  by  up  to  50,000 
times  while  maintaining  a 1:1  magnification  at  the  output  phosphor. 
The  two-stage  magnetically  focused  system  has  a 1000:1  dynamic 
range  and  high  signal-to-noise  ratio  (RCA  Electro-Optics  Handbook 
J974).  The  very  bright  output  image  permitted  the  use  of  commer- 
cially available  relay  lenses  for  focusing  the  image  with  1:1 
magnification  onto  the  photocathode  of  an  image  isocon  television 
pickup  tube  with  very  little  vignetting. 


The  image  isocon  (RCA  Electro-Optics  Handbook  3974)  video  camera  was 
selected  for  its  high  sensitivity  and  dynamic  range  (1000)  in  order 
to  match  the  dynamic  range  of  the  image  intensifier.  The  image 
isocon  operates  in  a 525  line,  60  f ield-per-socond  format. 

Additional  favorable  features  are  the  low  "lag"  (Ritman,  et  al. , 
1973) , low  noise,  ability  to  integrate  the  target  for  at  least  2 
seconds  (normally  1/60  second)  to  increase  the  "sensitivity",  to 
"zoom  in"  to  smaller  central  portions  of  the  image  and  to  tolerate 
very  bright  portions  within  the  image  (due  to  "raw"  beam  passing 
alongside  but  not  through  the  thorax)  without  loss  of  the  dynamic 
range  and  signal- to-- noise  within  the  image  of  the  thorax. 
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The  image  isocon  device  also  has  the  ability  to  retain  linearity  in 
regions  of  the  image  illuminated  at  x-ray  intensities  which  produce 
output  signals  below  the  "knee"  of  the  characteristic  brightness/sig- 
nal output  curve  while  simultaneously  a portion  of  the  x-ray  image 
intensity  exceeds  the  linearity  range  of  the  tube  so  that  a portion 
of  the  scene  is  imaged  beyond  the  normal  operating  range  of  the 
camera  (see  Fig.  6).  This  capability  is  particularly  important 
because  the  full  1000:1  dynamic  range  of  the  imaging  chain  can  be 
utilized  for  imaging  the  anatomic  structure  of  interest  (e.g.,  heart), 
when  portions  of  the  image  generated  by  unat. tenanted  "raw"  x-ray  beam 
image  are  intentionally  allowed  to  exceed  the  "knee"  of  the  char- 
acteristic: curve  of  the  imaging  chain.  This  practice  permits 
discrimination  of  x-ray  intensity  differences  down  to  four  orders 
of  magnitude  of  the  unattonuated  "raw"  beam  and  thereby  correspond- 
ingly increases  the  potential  density  resolution  of  the  cross- 


sectional  image  of  the;  structure  of  interest. 


The  reconstructed 


-images  obtained  from  this  type  of  scanning  are  subject  to  some 
distortion  confined  majorly  to  the  region  of  structural  border 
at  which  very  large  changes  in  x-ray  density  occur  (e.g.,  the 
ternal  surfaces  of  the  thorax)  and  which  have  no  apparent 
significant  efiect  on  the  quality  of  the  reconstruction  of  in 
ternal  structures  (e.g.,  the  heart).  Successful  methods  for 
correction  of  this  type  of  distortion  have  been  implemented 
(Ruesgegger  and  Ritman,  unpublished  data) . 
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Normally  the  video  image  is  recorded  on  video  tape,  stop-action 
video  disc  or  is  converted  directly  from  analog-to-digital  form 
for  storage  in  the  computer  (Robb,  et  al.,  1973). 

The  efficiency  of  the  imaging  chain  is  depicted  by  tracing  a single 

x-ray  photon  through  the  system  as  shown  in  Fig.  2.  The  absorption 

of  the  fluorescent  screen  is  approximately  50';  (Sturm  and  Morgan 

1 9 d S ) so  that  if  two  x-ray  photons  are  incident  on  the  screen  only 

one  will,  on  the  average,  be  absorbed.  If  the  x-ray  continuous 

spectrum,  i.e.,  the  bremsstrahlung  energy,  is  integrated  and  the 

photons  are  normalized  to  70  kev  photons,  at  normal  fluoroscopic 

levels  and  x-ray  to  fluorescent  screen  distances,  there  would  be 

7 2 

approximately  2 x 10  photons  per  mm  ' incident  on  the  anatomic 
structure  per  second. 


One  of  these  70  kev  photons  absorbed  by  the  screen  would  produce 
approximately  5,000  light  photons  (Sturm  and  Morgan  1949) , of  whit n 
501  would  be  useful  on  the  camera  side  of  the  screen.  The  efficiency 

I 

of  the  specially  coated  objective  lens,  50  mm  focal  length  f/0.75  is 
0.2%,  therefore  approximately  13.5  light  photons  would  be  projected 

V j 

onto  the  input  photoelectric  surface  of  the  image-intensif ier . Thu 
efficiency  of  this  surface  is  11%  - so  that,  this  one  x-ray  photon 

Ik  I 

* 


absorbed  by  the  screen  produces  1.4  photo  electrons.  Therefore, 
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this  becomes  the  weakest  link  in  the  chain  since  from  this  point 
on,  the  photoelectrons  are  accelerated  to  impinge  on  a phosphor 
producing  a large  increase  in  light  photons.  In  fact,  this  in- 
tensifier,  together  v.’ith  the  image-isocon , is  normally  used  with  a 
gain  of  more  than  25,000  in  equivalent  light  intensity. 

What  does  this  mean  in  terms  of  television  system  performance?  Since 
each  picture  is  viewed  by  the  television  camera  target  at  60  times 

7 

(i.e.,  fields)  per  second,  the  x-rays  in  the  raw  beam  (i.e.,  2 x 10 

2 5 

per  second  per  mm  ) m this  unit  is  reduced  to  about  3.3  x 10  x-ray 

2 

photon  per  field  per  mm  due  to  this  short  exposure  time  and  the 
increased  x-ray-to-subject  distance  as  shown  in  Column  A of  Fig.  2. 
Since  some  regions  of  the  thorax  (i.e.,  the  heart,  spine,  and  chest 

walls)  absorb  practically  992  of  this  radiation,  this  reduces  to 

3 2 

3.3  x 10  photons  per  field  or  exposure  per  mm'  shown  in  Column  13. 

The  number  of  photons  imaged  per  unit  area  can  be  increased  if  the 
minimal  area  of  interest  is  increased.  Thus,  by  averaging  the  image 
signal  from  4 adjacent  video  lines  (lilies  about  1 mm  apart)  within 
the  image  approximately  a two-fold  increase  of  signal- to-noise 
(i.e.,  A)  can  bo  achieved.  The  validity  of  this  approach  depends 
on  the  "independence"  (i.e.,  statistical  nature)  of  the  noise  in 
the  image  chain  output  signal. 
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The  statistical  nature  of  the  imaging  chain  noise  has  been  verified. 
Averaging  of  signal  data  samples,  either  taken  at  successive  points 
in  time  from  the  same  point  in  the  image,  or  from  contiguous  regions 
in  a uniform  image  obtained  at  one  point  in  time,  increases  the 
signal-to-noisc  in  proportion  to  the  square  root  of  the  number  of 
averaged  data  points.  Fig.  3 illustrates  this  finding. 

The  ability  to  "drive"  the  raw,  unattenuated  x-ray  beam  portion 
of  the  image  over  the  "knee"  of  the  characteristic  curve  of  the 
image  isocon  camera  can  be  utilized  to  increase  the  signal  in 
that  section  of  the  image  that  is  of  interest  (e.g.,  the  heart) 
by  integrating  images  generated  by  sequential  x-ray  pulses  on  the 
camera  tube  target.  This  is  accomplished  by  disabling  the  60-per- 
second  electronic  sweep-off  mechanism  until  an  adequate  signal  has 
been  accumulated.  Fig.  A illustrates  how  the  video  tube  target 
integration  time  can  be  controlled  via  the  computer.  This  approach 
effectively  increases  the  radiation  flux-per-image  pixel  while 
retaining  normal  x-ray  radiation  flux  per  pulse.  This  procedure 
reduces  the  number  of  images-per-second  that  can  be  recorded  in 
direct  proportion  to  the  number  of  successive  60-per-second  x-ray 
images  integrated  on  the  target.  This  technique  has  also  been  used 
to  image  the  distribution  of  radioisotopes  on  the  basis  of  their 
emitted  radiation  (Berggren,  ct  al.,  1976).  Up  to  a ten-fold 
increase  (relative  to  the  radiation  required  in  a small  chest)  in 


-9- 


x-  or  ganxna-ray  f lux-per-image  pixel  may  be  necessary  to  adequately  image 
the  heart  in  a large  chest.  This  will  require  up  to  150  msec  "inte- 
gration" times  between  successive  "read  outs"  of  the  isocon  camera 
target . 

ig . 4 

RESULTS 

a ) Operational  Charactcristi  cs 

Due  to  the  use  of  a flat  fluorescent  screen  and  flat  photo- 
cathodes of  the  image  .intensifier  and  image  isocon  tubes,  very 
little  geometric  distortion  occurs  over  the  entire  field  of  view 
as  shown  in  Pig.  5a.  The  lenses  and  image  isocon  systems, 
however,  exhibit  radial  intensity  fall-off  as  shown  in  the 
intensity  plot  of  a single  selected  video  line  (brightened  in 
image)  as  shown  in  panel  c.  The  image  brightness  at  the  peri- 
phery is  00%  of  the  brightness  at  the  center.  The  dynamic 
range  and  signal- to-noise  characteristics  are  illustrated  in 
panels  b and  d. 

ig.  5 

The  system  input/output  linearity  is  demonstrated  in  Fig.  6.  In 
the  presence  of  monochromatic  x-ray  and  absence  of  imaged  x-ray 
scatter,  the  Beer-Lambert  law  should  be  obeyed  (Ileintzcn  and 
Moldenhauor  1971)  so  that  the  logarithm  of  the  image  intensity 
should  linearly  decrease  with  increasing  depth  of  water  tra- 
versed by  the  x-ray  beam. 


g.  6 
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Although  operation  of  the  imaging  chain  so  that  segments  of 
the  image  exceed  "the  knee  of  the  curve"  as  illustrated  in 
curve  II  results  in  loss  of  some  information  (e.g.,  the  outer 
chest  wall)  the  resultant  increase  in  ability  to  discriminate 
small  changes  in  x-ray  intensity  within  the  regions  of  greatest 
interest  makes  improved  roentgen  density  resolution  within 
the  cross-sectional  images  of  these  regions  (e.g.,  the  heart) 
possible.  That  operation  of  the  image  isocon  video  camera  in 
the  integrating  mode  provides  an  increase  of  the  image  pixel 
signal  which  relates  directly  to  the  number  of  pulsed  x-ray 
images  integrated  is  shown  in  Fig.  7.  This  linear  relationship 
is  maintained  in  the  presence  of  "saturation"  if  the  image 
signal  in  certain  regions,  such  as  when  the  raw  beam  is  imaged 
alongside  the  chest. 


b)  Application  of  the  System 


A conventional  projection  image  of  a canine  thorax  and  one 
selected  (out  of  a possible  250)  cross-sectional  image  computed 
from  the  projection  images  of  the  chest  at  28  angles  of  view 
around  100°  (Robb,  et  al.,  1976)  is  shown  in  Fig.  0.  On  the 
left,  the  raw  beam  outside  the  chest  and  the  catheter  inside 
the  heart  are  clearly  detected.  Pulmonary  vasculature  and  the 
inferior  vena  cava  are  clearly  visible  within  the  chest.  The 
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limited  number  of  angles  of  view,  the  effect  of  x-ray  scatter 
and  the  mechanical  imprecision  of  the  rotation  of  the  dog  in 
front  of  the  imaging  chain  .results  in  a true  spatial  and 
density  resolutions  of  6 mm  and  105,  respectively,  within  the 
cross  section.  Contrast  medium-filled  coronary  arteries  of  less 
than  1 mm  diameter  can  be  visualized  in  the  intact  heart  when 
the  image  isocon  is  operated  so  that  those  regions  of  the  image 
with  low  information  content  (raw  beam  and  chest  wall)  are  over 
the  "knee"  of  the  characteristic  brightness/signal  output  curve 
of  the  camera,  scattered  radiation  is  reduced  by  collimination 
of  the  incident  and  transmitted  x-rays  to  a fan  beam  and  a 
focused  colliminator  oriented  at  right  angles  to  the  plane  of 
this  fan  beam  is  used. 


DISCUSSION 

| 

This  prototype  x-ray  imaging  system  has  been  shown  to  be  adequate 
for  three-dimensional  cross-sectional  reconstruction  of  the  intact 
thorax  and  its  contents.  The  anatomic  shape  and  dimensions  of  the 
lungs  and  epicardial  surfaces  can  be  imaged  without  the  use  of 
intravascular  agents.  These  GO-per-sccond  images  can  be  viewed  in 
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the  epi-  and  endocardial  surfaces  of  the  heart  at  60  images/second 
can  be  achieved,  if  at  least  101  concentration  of  iodinated  roentgen 
contrast  medium  is  present  in  the  blood  stream.  Reconstruction  of 
the  epicardial  and  endocardial  surface  will  permit  quantitative 
analysis  of  cardiac  anatomy  and  function.  In  addition,  the  presence 
of  501  concentration  contrast  medium  in  smaller  vessels  such  as  the 
coronary  arteries  will  permit  quantitative  analysis  the  three-dimen- 
sional anatomy  and  segmented  flow  through  the  coronary  macrovasculature 
and  regional  myocardial  perfusion.  Detection  of  disease  processes 
which  perturb  the  anatomic  structure  or  function  of  the  circulatory 
system  in  any  region  of  the  body  will  also  be  possible. 

Use  of  this  prototype  imaging  chain  in  a reconstruction  (computerized 
tomography)  system  is  a compromise  between  speed  of  data  acquisition 
(i.e.,  temporal  resolution)  and  spatial-density  resolution.  Future 
modifications  of  the  system  should  be  aimed  at  increasing  the 
efficiency  of  the  imaging  chain  and  reduction  of  degradation  of  the 
image  due  to  x-ray  scatter  and  other  efforts  to  increase  the  signal- 
to-noise  ratio  of  the  overall  system.  Improvement  of  the  x-ray 
focused  grid  and  fluorescent  screen  are  especially  required. 
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LEGEEDS 

Fig . 1 Schematic  arrangement  of  a single  roentgen  video  imaging 
chain  for  high  temporal  resolution  computerized  tomography,  a dynamic 
spatial  reconstruction  system.  The  intensity  of  each  x-ray  pulse, 
every  1/60  second,  is  monitored  by  the  ionization  chamber  so  as  to 
provide  a value  proportional  to  radiation  incident  on  the  subject- 
under  study  (see  text  for  details).  (Cardiovascular  Imaging  and 
Processing.  Theory  and  Practice.  Soc.  Photo-Optical  Inst.  Eng. 
72:103-122  (1975)) 

Fig.  2 Efficiency  of  imaging  chain  depicted  by  tracing  x-ray  photons 
through  system.  A theoretical  relationship  of  the  signal-to-noise 
of  this  type  of  imaging  chain  has  been  derived  and  is  expressed  by 
the  equation  for  S/N.  This  equation  gives  a theoretical  signal-to- 
noise  ratio  of  35/1  in  the  thickest  part  of  the  chest  (i.e.,  where 
greatest  x-ray  absorption  occurs)  (see  text  for  details).  (Pro- 
ceedings of  the  San  Diego  Biomedical  Symposium,  February  4-6,  1976, 
San  Diego,  California,  in  press) 

Fig.  3 Signal  noise  in  the  SSDSR  imaging  chain  was  measured  as 
four  standard  deviations  of  the  samplo-to-sample  variation  of  the 
voltage  of  the  video  signal  about  the  mean  value  measured  at  the 
same  site  (i.e.,  pixel)  on  a single  video  line.  Panel  (a)  - the 


intensity  of  the  highly  collimated  x-ray  beam  was  attenuated  by 
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placing  10  and  20  cm  water  between  the  x-ray  source  and  fluorescent 
screen  while  KV  and  ma  of  the  x-ray  system  were  left  unchanged. 

The  dashed  curve  was  calculated  on  the  theoretical  basis  that 
the  noise  in  the  video  signal  should  be  proportional  the  square 
root  of  the  magnitude  of  this  signal.  The  magnitude  of  the  cal- 
culated noise  was  assumed  to  be  identical  to  the  measured  noise  at 
the  700  mv  signal  level.  The  measured  noise  falls  less  rapidly 
than  the  calculated  values  indicating  that  noise,  over  and  above 
that  to  be  expected  from  the  x-ray  image,  contributes  to  the  imaging 
chain  signal  noise. 

Panel  (b)  - The  measured  signal  noise  (solid  lines)  decreased 
similarly  to  the  calculated  (dashed  lines)  when  an  increasing 
number  of  signal  samples  were  averaged.  The  calculated  and  mea- 
sured values  are  for  a video  signal  level  of  700  mv  measured  at  the 
same  site  on  the  same  video  lino  (i.e.,  the  same  pixel)  at  successive 
30/sec  instants  in  time  (curve  I) . The  measured  decrease  in  noise 
is  virtually  identical  to  the  calculated  values,  indicating  that 
the  noise  is  statistical  in  nature.  When  the  video  signal  values 
were  measured  at  the  same  site  on  different  successive  video  lines, 
i.e.,  from  the  same  vertical  column  of  A/D  samples,  in  one  image, 
the  statistical  nature  of  the  noise  was  also  confirmed  (curve  II). 
Curve  III  is  based  on  values  obtained  by  averaging  signal  samples 
along  the  same  horizontal  video  line  in  the  same  uniform  gray 
level  image.  This  curve  indicates  that  there  is  a non-statistical 
bias  in  the  noise  (possibly  caused  by  variations  in  the  structure 
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bias  in  the  noise  (possibly  caused  by  variations  in  the  structure 
of  the  image  or  signal  induced  oscillations  in  the  electronic 
circuitry. 

Fig.  4 Schematic  of  assembly  to  provide  programmable  computer 
control  of  the  sweep  (read-off)  frequency  of  a video  camera.  Be- 
cause of  the  integrating  characteristics  of  the  input  imaging 
target  of  the  camera,  when  high  temporal  resolution  is  not  required, 
its  sensitivity  to  low  light  levels  can  be  increased  by  decreasing 
the  sweep  frequency  from  the  normal  30  or  60-per-second  value  to 
1-pcr-second  or  less  (see  Fig.  7) . 

Fig.  5 a.  Photo  of  x-ray  image  of  test  grid  consisting  of  crossed 
radiopaque  wires  2 cm  apart,  and  cross  hatchings  2 mm  apart  along 
the  central  cross  wires.  The  slight  barrel  distortion  in  this 
photograph  is  primarily  due  to  the  videomonitor  because  the 
horizontal  wires  remain  in  alignment  with  the  horizontal  video 
lines. 

b.  Photo  of  x-ray  image  obtained  of  a step-wedge  containing  1,  8 
15,  22,  and  29  layers  (each  0.63  cm  thick)  of  pressed  wood  using 
120  kv  peak  x-ray  beam  pulsed  60-per-second,  each  pulse  0.34  msec 
duration.  Two  steel  wires  0.5  and  1.0  mm  diameter  were  strung 
diagonally  across  the  wedge.  In  this  system,  107,  of  the  incident 


x-ray  beam  was  transmitted  through  18  layers  of  pressed  wood  as 
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measured  at  the  fluorescent  screen  with  a Victoreen  radiation  meter. 
The  video  signal  deflection  (panel  d)  caused  by  the  1 mm  steel  wire 
is  clearly  distinguishable  in  the  video  signal  corresponding  to  the 
22  layers  of  pressed  wood.  Direct  observation  of  the  x-ray  image 
on  the  videomonitor  permits  clear  visualization  of  the  1.0  mm  wire 
through  the  29  layers  of  pressed  wood. 

c.  The  video  signal  of  a single  horizontal  video  line  near  the 
middle  of  the  roentgen  video  image  in  panel  a is  displayed  to  il- 
lustrate the  radial  fall-off  in  image  brightness.  The  amplitude  of 
the  signal  above  zero  is  directly  proportional  to  image  brightness 
and  the  square  wave  deflections  on  the  far  left  and  right  are  video 
synchronization  signals.  Most  of  the  fall-off  is  attributable  to 
the  optical  system  and  ihe  operational  characteristics  of  the  image 
isocon. 

d.  The  video  signal  of  a single  horizontal  video  line  near  the 
middle  of  the  roentgen  video  image  in  pane]  b is  displayed  to 
illustrate  the  decrease  in  video  signal  noise  (thickness  of  trace) 
with  decrease  in  signal  magnitude.  The  peak-to-pcak  noise  in  the 
video  image  is  approximately  10%  of  peak  signal  ("raw"  beam)  and  1% 
over  the  29  layers  of  pressed  wood  where  the  signal  is  about  1%  of 
"raw"  beam.  The  deflections  caused  by  the  steel  wires  are  visible 
in  the  signal  corresponding  to  the  x-ray  transmitted  through  the  22 
and  15  layers  of  pressed  wood  for  the  1 mm  and  0.5  mm  steel  wires, 
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respectively . (Proceedings  of  the  San  Da  ego  Biomedical  Symposium, 
February  4-6,  1976,  San  Diego,  California) 

Fig.  6 Increasing  depth  of  water  transradiated  by  100  kv  x-ray  beam 
results  in  exponential  decrease  in  video  signal  as  shown  in  curve  I. 
At  increasing  depths,  the  decrease  of  video  signal  becomes  less  rapid 
due  to  "hardening"  of  the  remaining  x-ray  beam.  The  relationships  in 
curve  II  were  obtained  by  increasing  the  x-ray  beam  intensity  so  that 
the  signal  from  the  x-ray  beam  attenuated  by  less  than  5 cm  of  water 
exceeded  the  linear  segment  (the  "knee")  of  response  curve  of  the 
image  isocon  camera. 


Fig.  7 Integration  of  successive  pulsed  x-ray  images  on  the  image 
isocon  tube  target  results  in  a linear  increase  in  the  video  signal 
amplitude  when  this  is  "read"  out  after  completion  of  the  integration 
period.  To  prevent  non-linear  integration  at  very  low  image  in- 
tensities, the  periphery  of  the  x-ray  image  was  highly  illuminated 
during  the  integration  period  so  as  to  maintain  proper  target 
potential.  The  noise  in  the  imaging  system  remains  proportional  to 
the  video  signal  (see  Figure  3a)  so  that  a ten-fold  increase  in 
signal  from  70  mv  to  700  mv  results  in  an  increase  in  noise  of 
approximately  six-fold  (theoretically  approximately  three-fold)  so 
that  a two-fold  increase  of  signal  noise  results  in  addition  t.o  the 
ten-fold  (more  than  3 bit)  increase  in  analog- to-digital  conversion 
resolution  in  that  part  of  the  image  of  interest. 
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Fig.  8 Loft  panel:  Monoplane  lateral  roentgen  projection  image 

of  dog  thorax.  The  brightened  line  indicates  the  level  of  the 
reconstructed  cross  section  shown  in  right  panel. 

Right  panel:  Cross-sectional  image  was  computed  from  roentgen 

density  profiles  measured  from  same  horizontal  video  line  recorded 
at  28  equispaced  angles  of  view  around  an  arc  of  160°.  Similar 
cross-sectional  images  can  be  computed  from  each  of  the  240  hori- 
zontal video  lines  contained  in  each  60-per-second  video  field 
image.  Physiologic  stationarity  of  the  respiratory  and  cardiac 
cycles  plus  synchronization  with  the  x-ray  pulsing  and  video  record- 
ing cycles  are  required  for  accurate  reconstruction  of  the  thorax 
of  intact  anesthetized  dogs  when  a single  roentgen-video  image 
chain  and  consequent  mandatory  mechanical  circumferential  scanning 
are  used  (R.  E.  Sturm,  et  al..  Proceedings  of  the  San  Diego  Bio- 
medical Symposium,  February  4-6,  1976,  San  Diego,  California,  in 


press) 
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